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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 

DAMPING IN PITCH OF LOW-ASPECT-RATIO WINGS 
AT SUBSONIC AND SUPERSONIC SPEEDS 
By Murray Tobak 


SUMMARY 


The concept of indicia! functions is applied to the analysis of the 
aerodynamic phenomena associated with the short-period pitching mode of 
wings in subsonic and supersonic flight. Simple physical relationships 
are pointed out and are used to study the effect on the rotary-damp ing- 
moment coefficient of changes in center-of -gravity position, Mach number, 
aspect ratio, plan form, frequency, and thickness. Qualitative conclu- 
sions are drawn from the results of this investigation and are compared 
with the results of experiments for a series of low-aspect-ratio wing- 
body combinations having triangular, swept, and unswept wing plan forms . 

Results of the experimental Investigation, which were obtained by a 
single-degree-of-freedom free-osci.] lation technique over the Mach number 
ranges 0.6 to 0.9 and 1.2 to 1 . 9 , were In good agreement with the results 
of theoretical computations. The predictions of ranges of supersonic 
Mach number and center-of -gravity positions over which dynamic instabil- 
ity may be expected, of the beneficial effect on the damping moment of a 
reduction in aspect ratio, and of only a small effect of thickness on the 
damping moment were confirmed by the experimental results. 

The occurrence at high subsonic Mach number's of small-amplitude 
self-sustained pitching oscillations is noted, and a hypothesis is 
advanced for its explanation. 


INTRODUCTION 


In the classical study of the longitudinal motion of an aircraft, 
it is usually found that the motion resulting from a small equilibrium- 
destroying disturbance consists of two modes: one, a lightly damped, low- 
frequency motion at essentially constant attitude, called the phugoid 
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oscillation; the other, a rotary-pitching and plunging oscillation of 
high frequency (relative to the phugoid frequency) called the short- 
period oscillation. The phugoid oscillation has generally been described 
as resulting from a slow interchange of potential and kinetic energy as 
-.the aircraft experiences periodic variations in airspeed and altitude. 

The character of the phugoid motion as influenced by airspeed, altitude, 
and aircraft geometry has been well understood for some time (see, e.g., 
ref. l) . The short-period motion, on the other hand, having in the past 
been found to be highly damped and of short duration, has been the cause 
of no concern. Its characteristics therefore have not been as fully 
investigated as those of the phugoid oscillation. With the advent of 
flight at speeds approaching and exceeding the speed of sound, however, 
the loss of rotary danping occurring in practically all aircraft at 
speeds near the sonic speed has caused renewed interest in the short- 
period pitching mode. Unlike the easily controlled phugoid oscillation, 
the deterioration of damping in the short-period mode is of serious con- 
cern to the pilot, since the period of the oscillation can be of the 
same order of magn i tude as the pilot's reaction time. The oscillation 
may therefore be difficult or even impossible for the pilot to control 
manually. Furthermore, the additional load imposed upon the airframe 
due to a rapid growth of the amplitude of a negatively damped oscillation 
makes possible the occurrence of structural failure. It is therefore of 
considerable interest to obtain an understanding of the nature of the 
short-period mode, parallel to that which has been gained of the phugoid 
mode. 


One means of viewing the aerodynamic phenomena occurring during the 
short-period oscillation from a fundamental standpoint is through appli- 
cation of the concept of indicial functions. In this approach, the var- 
iations with time of the aircraft angle of attack and angular velocity 
during the oscillation are replaced by a large number of small instanta- 
neous or step changes. The transient aerodynamic reactions to these step 
changes are termed "indicial functions," and have been calculated theo- 
retically for several classes of wings (refs. 2 to 6) . By suitable 
superposition of these results (refs. 7 to 9)> the aerodynamic forces and 
moments caused by the given maneuver can be studied. It will be the 
purpose of this report to make such a study for the simplified case of an 
aircraft performing single-degree-of-freedom rotary oscillations. For 
this maneuver, which corresponds to the short-period oscillation when the 
plunging velocity of the aircraft 1 b zero, the use of simple physical 
relationships associated with the Indicial function concept enables 
qualitative studies to be made of the separate effects on the aerodynamic 
forces and moments of changes in Mach number, aspect ratio, plan form, 
frequency, and thickness. Results of this investigation are then con>- 
pared with the results of experiments with a group of low-aspect-ratio 
wing -body combinations. The tests were conducted in the Ames 6- by 6-foot 
supersonic wind tunnel and were similar in technique to those reported In 
reference 10. 



MCA EM A52L0lia 


3 


A 

c L 

Qm 

C p 

I 

Mq 

R 

S 

Vo 

a o 

b 

c 0 

c 

e 

i 

k 

Ap 


q. 

<lo 

t 

"^a 


NOTATION 


2 

aspect ratio, b /S 

lift coefficient, lift/c^S 

pitching-moment coefficient, pitching moment /^Sc 

P-Po 

pressure coefficient, 

<lo 

moment of inertia, slug-ft 2 
free-stream Mach number, V 0 /ao 

Reynolds number, based on wing mean aerodynamic chord 
wing area, including portion enclosed by body, sq ft 
flight speed, ft/sec 

speed of sound in free stream, ft/s ec 

wing span, ft 

wing root chord, ft 

2 f b/2 .2 

wing mean aerodynamic chord, — J (local chord) dy 

o o 

base of natural logarithms 


reduced frequency parameter, wc/2V 0 

local loading at plan-form surface, 

pressure lower surface - pressure upper surface 

angular velocity due to pitching, radians/sec 
free-stream dynamic pressure, p Q V 0 2 , TP/sq ft 
time, sec 


time required following an instantaneous change in angle of attack 
or angular velocity for the transient lift or moment to attain 
steady state, sec 
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x,y, z rectangular coordinates 

x a.c. distance from leading edge of M.A.C. to aerodynamic center, ft 
x Q distance from, leading edge of M.A.C. to axis of rotation, ft 

AXq Xa. c . - x o 

a angle of attack of wing center line with respect to free-stream 

direction (sketch a) 

0 yiMo 2 -ii 

7 ratio of specific heat at constant pressure to that at constant 
volume 

& airfoil-thickness ratio, maximum thickness /chord 

0 angle of wing center line with respect to horizontal axis 
(sketch a) 

V acute angle between wing plane of symmetry and trailing edge 
(sketch p) 

p Q free-stream density, slugs/cu ft 

to angular frequency of oscillation, radians/sec 

o a distance traveled, measured in half M.A.C. lengths, in the time 
interval t a , 2V 0 t a /S 

When a, a, and q are used as subscripts, a nopdimensional deriv- 
ative is indicated, and this derivative is evaluated as the independent 
variable (a, a, or q) approaches zero. For example, 


C % = 


"I 



c mi = 

dC m 

_ da J 

a-> o 

s a c 



"c>C m “ 





d sL. 

L 2V °J 

q— >o 



ct — ^ o 


0 
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THEORY 

Application of Indicial Functions to the Aerodynamic 
Theory of Unsteady Flows 


One of the most useful tools in the study of unsteady flows is the 
concept of indicial aerodynamic functions which may he defined briefly 
as the aerodynamic response of the airfoil as a function of time to an 
instantaneous change in one of the conditions determining the aerody- 
namic properties of the airfoil in a steady flow. Theoretical aerody- 
namic indicial functions were first derived by Wagner (ref. 2) for the 
two-dimensional wing in incompressible flow. More recently, these 
results have been extended by Heaslet and Lomax to cover the compressible 
case for both subsonic and supersonic speeds (ref. 4) . In addition, 
theoretical indicial functions have now been obtained for both wide and 
slender triangular wings and rectangular wings, all for supersonic speeds 
(refs. 4 to 6) . 

The Indicial function derives its usefulness primarily through the 
ease with which it lends itself to the powerful and well-established 
methods of the operational calculus (refs. 7 "to 9) • With the use of these 
methods, the aerodynamic forces and moments caused by arbitrary motions 
of the airframe can be studied from a fundamental standpoint. Because of 
the wide range of applicability of this means of approach in unsteady flow 
analyses, a considerable portion of the succeeding discussibn is devoted 
to the fundamentals involved. 


Definition of coordinate system .- In the succeeding analysis the 
stability system of axes is used. The origin of the coordinate system 
is placed in the airfoil so that the y axis which Is perpendicular to 
the vertical plane of symmetry is coincident with the axis of rotation of 
the airfoil; the positive branch of the x axis is pointed in the direc- 
tion of flight; and the z axis lies in the vertical plane of symmetry, 
positive downward. The angle of attack a is measured as the angle 
between the chord plane of the airfoil and the xy plane, and is shown 
as positive in sketch (a) . The angle of pitch 9 is the angle between 
the chord plane of the airfoil 
and the horizontal plane (an 
arbitrary reference) and is also 
shown positive in sketch (a) . 

Forces are measured as positive 
upward, whereas pitching moments 
are positive when tending to 
increase the angle of pitch in 
the positive direction. When 
the airspeed V Q is constant 
which corresponds to the 



Sketch (a) 
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condition under study, the translatory and angular motions of the air- 
foil with respect to any system of coordinates are defined if the time 
histories of the angle of attack a and the angle of pitch 9 and 
their derivatives are known. For purposes of clarity, two different 
harmonic motions of the aircraft are shown, illustrating the difference 
between a flight path which involves a constant angle of attack and a 
varying angle of pitch (sketch (b)) and one which involves a constant 
angle of pitch and a varying angle of attack (sketch (c)). 



ang/t of attack * O angle of attack * a 

Sketch (b) Sketch (c) 


Now consider the case of a wing executing harmonic rotary oscilla- 
tions about the y axis while the origin of the coordinate system 
traverses a level path at constant velocity V Q . This case corresponds 
to that of a wind-tunnel model mounted to permit single-degree-of- 
freedom rotary oscillations, or to the short-period mode of an aircraft 
in flight when the plunging velocity of the center of gravity is zero. 
Here a and 9 are equal, so that the maneuver is defined by one vari- 
able, the time history of either a or 9. Let the angle of attack 
be a and the angular velocity be q (q = dS/dt = da,/dt). At any 
instant, the normal velocity at any point on the airfoil surface is com- 
posed of. two parts, one due to the instantaneous angle of attack aV 0 , 
the other due to the angular velocity at the same instant -qx (see 
sketch (d)). These are two of the instantaneous boundary conditions of 
the unsteady flow. 



2X 


NACA RM A52L04a 


7 


Solutions for the aerodynamic forces and moments which correspond 
to these boundary conditions may be derived by a number of methods 
involving various degrees of approximation. In succeeding sections, 
the use of the concept of indicial functions and the principle of super- 
position for this purpose will be illustrated and compared with other 
current widely used methods. 

Concept of indicial functions .- In order to illu strate this concept, 
assume that the airfoil under consideration has been flying a level path 
at zero angle of attack. At some time, which is designated time zero, 
the wing is caused to attain simultaneously a constant angle of attack a 
and angular velocity q. The normal velocity of the flow next to the 
surface of the airfoil therefore changes discontinuously from zero to a 
pattern that is constant with time and identical in shape to the pattern 
shown previously in sketch (d) . The lift and pitching moment that 
result are of a transient character and attain their steady-state values 
corresponding to these new boundary conditions only after a significant 
interval of time has passed. It should be noted there exists an essen- 
tial difference between the length of this time interval at subsonic and 
supersonic speeds. At supersonic speeds, the vorticity shed into the 
airfoil wake cannot influence the flow about the airfoil but at subsonic 
speeds this influence exists for all time. The result is that the lift 
and moment reach steady-state values in a finite time at supersonic 
speeds but approach these values asymptotically at subsonic speeds. In 
either case, however, the time responses in lift and moment to the step 
changes in a and q. are termed indicial functions. Sketch (e) illus- 
trates typical subsonic and supersonic indicial lift responses to a step 
change in the angle of attack. 




Sketch (e) 
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It is obvious that tlie time history of the wing motion during a _ 
short-period oscillation may be broken down into an infinite number of 
infinitesimally small step changes in the angle of attack and step 
changes in the angular velocity. The summation of the indicial lift and 
moment for these steps then yields the total lift and moment at any pre- 
scribed time. In sketch (f), the mechanics of the procedure are illus- 
trated for an arbitrary angle-of -attack variation. Here, the given 

angle-of -attack variation Ib 
replaced by a number of small 
step changes. Within each step 
the corresponding response in 
lift is shown plotted for con- 
venience. It is then apparent 
that the total lift at time t 
is equal to the sum of the 
increments of lift in each step 
at time t. As indicated by 
the leaders, however, it is 
clear that the increments of 
lift for the various steps at 
time t are equivalent to 
increments in the first step 
at time t - t 1 . Alterna- 
tively, then, the total lift 
at time t can be written as 1 


Sketch (f) 

t 

Cl(*) = CL a (t)a(0) Cx^Ct-ti) || (t 1 )At 1 (l) 

o 

After a transformation of variables, t - t^ = T , and letting the incre- 
ment of time approach zero, equation (l) canjbe rewritten.in a form of 
Duhamel ' s integral (see, e.g., ref. 9) 

C L (t) = A J* Cl jT)a(wr)dT (2) 

1 Here, and in the remainder of this report, the use of parentheses Is 
reserved solely for the indication of functional dependence. Thus, 
for example, in equation (l) the term Cj^Ct-ti) Is interpreted as 
the value of Cj^ at time t - t x , whereas ~ (t x ) is the value 

of Ax /At at time t All other enclosures indicate algebraic 
expressions in the usual sense. " 
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A similar procedure is carried out for the angular velocity 
variation, whereupon the total lift coefficient at the prescribed 
time t becomes 

C L (t) ' -k C C L a < T > a < t - T > aT + 4 -C %< T > 2^ l(t ' T)dT <3) 

It should be pointed out that in this form equation ( 3 ) is applicable 
to the analysis of arbitrary motions, the only restriction being that 
the flight speed is constant. In the following sections, however, the 
application of equation ( 3 ) is restricted to harmonic motions having a 
single degree of freedom. The reasons for this restriction are two- 
fold: first, the motions of a statically stable aircraft in response to 
a disturbance are most generally of a harmonic nature; and second, such 
a restriction permits an assessment of the influence of the time rate 
of the airfoil motions on the aerodynamic forces and moments. 

Application of Indicial functions to harmonic pitching oscillations 
Consider first a pure sinusoidal pitching oscillation, the angle of 
attack being zero throughout the motion. The flight path for such a 
motion has been illustrated in sketch (b) . In this case, the angle of 
pitch is given by 

e(t) = 0 o ^ ia)t 

where 0 O is the maximum amplitude of oscillation and w is the 
angular frequency. The angular velocity is, of course, q(t) = 0 = 
ito0 o e iajt = ico0(t) . Inserting the value for q(t) in equation ( 3 ) and 
performing the indicated operations, there results 


°lW - - e(t) 4 %<T) e -i» T dT + Soc^t) 


1U)C 


(*0 
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Note in sketch (g) that Ci^(t) is equal to Cj^^t) -F 2 (t) , and that for 
subsonic speeds F 2 (t) approaches zero as T approaches t. 



Sketch (g) 

Replacing CLq(t) in equation (4) by this equality. 


= h5 o, (t) + £i ; t F 2 (T)e 
e(t) 2V 0 ^ 2V 0 J o 2 


-iUT 


dr 


(5) 


For subsonic speeds, let t approach infinity. With this sub- 
stitution, equation ( 5) thereby represents the lift coefficient due to 
the harmonic pitching motion after the transient loading subsequent to 
the start of the motion has reached steady state. Then separating 
equation ( 5 ) into components in-phase (real part) and out-of-phase 
(imaginary part) with 0, there is obtained 


— = f F 2 (t) cos cutIt 
0 2V 0 j o 2 

-t = Cl (») - u / Fo(t) sin WTdT 

iwc0/2V o V J o , 


( 6 ) 


Introduce the nondimens ional parameters. 


<P 



number of half M.A.C. lengths traveled in time t 


k 


cue 

2V^ 


reduced frequency 


y 


-O 



mCA RM A52L04a 


11 


In terms of these parameters, equation ( 6 ) becomes, for Mq<1, 

— = k 2 * * * f*° F 2 (cp) cos kcpdcp 
0 J o 

— 7 — = ClqCc) - k /“ F a (q>) sin kcpdcp 
qc/2V 0 ^ o 


(T) 


At supersonic speeds, equations (7) may be simplified somewhat 
since the build-up in lift is completed in a finite number of half 
M.A.C. lengths of travel cr a . In equations (l) , therefore, the upper 
limits of the integrals may be replaced by cr a , since beyond that point 
F 2 (cp) is zero. 


— = k 2 f a F 2 (<p) cos kcpdcp 

9 ° cp > 0 , 

M 0 > 1 

= C Lq (a a ) - k / & F 2 (q>) sin kcpdcp 


( 8 ) 


Thus, it appears from equations (7) and ( 8 ) that there are both in-phase 
and out-of-phase lift forces associated with the harmonic pitching oscil- 
lation. Notice, however, in equations ( 8 ) that if the trigonometric 
terms are expanded and the reduced frequency is required to be small 
compared to unity (corresponding to the frequencies encountered in 
dynamic stability work) terms containing second and higher powers of k 
will be very small compared to first-order terms. Thus, for slow fre- 
quencies, the only force of consequence during the pitching oscillation 
is the first order in frequency out-of -phase lift force, qc/2V 0 C]^(a a ). 2 


2 This quantity is, of course, the same lift force in phase with the 

pitching velocity which would occur alone had the wing been executing 

a steady turn (q constant). One of the chief advantages of the indi- 
cial response method, at least for supersonic speeds, is the ease with 

which the relative importance of the various terms contributing to the 

total lift and moment can be assessed and the sources of the important 
contributions identified. 
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qSC L (cr g )^ 

2V 0 


CJlota!) 


~ Si k k<*)sin Wm 9 k*/i*(4)cosk$d$ 

2V 0 Jo *4 


Sketch (h) 


The phase relationships for the har- 
monic pitching oscillation are indi- 
cated in sketch (h) . It is evident 
from sketch (h) that the total lift 
leads the angle of pitch by nearly 90 ° 

Application of indicial functions 
to harmonic plunging oscillations . - 
Next, consider a purely sinusoidal 
variation of the angle of attack, the 
angle of pitch being zero throughout 
the motion. The flight path for this 
motion has been illustrated in 
sketch (c). Here, a, equals ccoe iw , 
where, as previously, ao amd. w are 
the maximum amplitude and angular fre- 
quency, respectively. Applying equa- 
tion (3) again. 


C L (t) = ao ± £ CLa (T)e iU)[t - T] dT 


= itoa(t) CLa(T)e” iuT dT + ^oCx^t) (9) 

Now, as in the previous example, let C^(t) = Cj^Ct) - Fi(t) so that 
equation (9) becomes 


C L (t) = a(t)Cn (t) - iwa(t) / Fi(t ) e" iWT dT 

o 


( 10 ) 


Again, introduce the nondimensional parameters 9 and k, separate 
equation (10) into its real and imaginary parts, and let 9 approach 
infinity for subsonic speeds and a a for supersonic speeds. There 
results 


TT ' c Iu< x > ' k 


sin kcpdcp 


Cl 

ac/2V 0 


f F x (9) cos k939 


Mo<li 
X = o a , M 0 >1> 


( 11 ) 


J 
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Notice in equations (ll) for supersonic speeds that when the trigono- 
metric terms are expanded for the slow frequency case, as was done in 
the previous example, there appears an in-phase term of zero order 
in k, aCi^Caa), and an out-of -phase term of first order in k. 


- _2£L f “ F 1 (cp ) dcp . These, then, are the principal contributions to the 

lift forces for the slow- frequency angle-of -attack variation. The phase 
relationships for this motion are shown graphically in sketch (i). For 
this case, it is evident that the total-lift force can lag behind the 
angle of attack because of 
the negative out-of-phase 
contribution. 


,°a 


ac r 

2V n 


F x (q)) cos kfpdcp 



Application of indicia! 
functions to harmonic rotary 
oscillations .- Finally, con- 
sider the case of harmonic 
rotary oscillations. Here, 
as previously mentioned, the 
normal velocity over the 
wing surface is composed of 
contributions from both the 
angular velocity and the 
instantaneous angle of 
attack, so that the complete 

expression in equation ( 3 ) must be employed to obtain the total lift. 
However, for single-degree-of -freedom rotary oscillations, a equals 9 
and a equals q, so that in this case the separate expressions given 
for the harmonically pitching wing (eqs. ( 7 ) and (8)) and the harmoni- 
cally plunging wing (eq. (ll)) can be combined to give the total lift 
for a wing executing harmonic rotary oscillations. Then, adding the 
results of equations (7), (8), and (ll), the in-phase and out-of -phase 
components of the total lift become 


w - c ^ (x) 


k F,(<p) sin kfpd9 + k 2 /^ F 2 (<p) cos k<pd9 
o o 


(12a) 


A=», Mo<l 
A.=o a , Mq> 1 


7 ~ — = C Ln (X) - k / F 2 (<p) sin kCpd9 - / X F,(<P) cos 

qc/2V 0 ^ 0 ° 


k9d9 


(12b) 
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The phase relationships for the rotary oscillation may, of course, also 
be obtained by directly adding the results given in sketches (h) and (i). 
The result of this addition is shown in sketch (j). 



It will be noted in sketch (j) that the total-lift force can either lag 
behind or lead the angle of attack, depending on the relative magnitudes 
of the three terms comprising the out-of-phase lift. The total lift is 
shown lagging behind the angle of attack in sketch (j), which situation, 
for axis positions ahead of the point of concentration of the total lift, 
gives rise to the possibility of the development of negatively damped 
rotary oscillations. 

Again, the complete frequency-dependent equations for the total 
lift of a wing in supersonic flow due to the rotary oscillation 
(eqs. (12), \ = cr a ) may be reduced to first order in k for the slow- 
frequency case in the same manner as was described in the two previous 
examples to give 


~ - < -'L a ( °a) 


CL 

qc/2Vo 


CL a (o a ) - / ff& F x (cp)dcp 
^ o 


( 13 ) 



NACA EM A52L04a 


15 


For all three examples, the same procedure may, of course, he used 
to obtain the pitching -moment coefficient. Only the pitching-moment 
equations for the rotary oscillation are presented below, since the 
correspondence between the lift and moment equations is obvious. For 
the rotary oscillation case, then, 


~ = C~(\) - k / X F S (<P) sin k<Pdcp + k 2 f X F 4 (<P) 
a “a* *o o 

q -y| j~ = Cn^CX) -k f* F 4 (q>) sin kspdcp - f X F s (<p) 


cos k<Pd9 1 


X=», 

X=Og_, 


Mq< 1 
Mq> 1 


(d 


cos kspdcp 


where, as previously. 


F s (<p) = Cmad ‘ 

and 

= Cj^Cx) - 0^(9) 


Again, reducing equations (l4) for supersonic speeds (X=o a ) to first order 
in frequency, there results 


TT = W ff a) 

> 


Cm 

qc/2V 0 


= ^(ffa) - f o ° a F s (cp) dcp 


(15) 


The complete frequency-dependent equations for the lift and pitching- 
moment coefficients for the rotary-oscillation case (eqs. (12) and (l4) ) 
describe completely the aerodynamic forces and moments resulting from 
the single-degree-of-freedom pitching mode. For the purposes of the 
present discussion, however, it is sufficient to limit consideration to 
the simpler first order in frequency results of equations (13) and ( 15 ) • 
The si gni fi can ce of the effect of the higher-order terms on the out-of- 
phase pitching -moment will be examined in a later section of this report. 


With regard to the first-order results for the supersonic lift and 
pitching -moment coefficients, it is instructive to note that the 
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quantities £ F 1 (q>)dcp and £ F a (cp)d<p in equations (13) and ( 15 ) 

are represented geometrically by the areas of the shaded portions of 
sketch (k) . 




Sketch (k) 

The manner in which these areas are affected by variations in Mach 
number, aspect ratio, plan- form shape, and thickness will be used as a 
guide in later sections of this report to determine the significance of 
these parameters. 

In the foregoing discussion, no mention has been made of reducing 
the complete equations for the lift and moment coefficients at subsonic 
speed (eqs. ( 12 ) and (l4), A,=oo) to first order in frequency as was done 
for the equations noted as applying at supersonic speed. It is evident 
that if the same procedure had been applied for subsonic speeds, the 
area corresponding to the term f F 1 (9)<1 C P can either be finite or can 
become infinitely large, depending on the manner in which the indicial 
lift function Cj^vCp) approaches its steady-state asymptote as cp->oo. 

In the latter case, there exists the interesting anamoly of an infinite 
out-of -phase lift force as the frequency approaches zero. As can be 
seen from the results of reference 4, such will be the case for the two- 
dimensional. wing. This result as the frequency approaches zero is not 
peculiar to the indicial analysis alone, but has been pointed out by a 
number of authors using different approaches. As indicated by Miles 
in reference 11 , however, the anamolous result can be considered to be 
a consequence of assuming a two-dimensional flow, and there is reason 
to believe that the difficulty as. the frequency approaches zero will not 
exist for finite-span wings. 

As has been mentioned previously, the use of the characteristic 
areas J a F 1 (cp)<i c P / a F 3 (tp)dcp will be shown to be of considerable 

value in estimating the <3amping-in-piteh characteristics of wings at 
supersonic speeds. For the two-dimensional wing at subsonic speeds, the 
singularity as k -^0 prevents the use of such a simplified approach 

GSBFmBHTMB* 
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without further study. However, rather than return to the use of the 
full frequency-dependent equations ((12) and (14)), the reduction of 
the equations for the out-of -phase lift and moment to first order in 
frequency will he made in such a manner as to preserve the significance 
of these areas. To accomplish this end, equation (l2b) is reconsidered. 
It is evident that the first integral in equation (12b) may be discarded, 
since its contribution to the out-of -phase lift is at least of second 
order in frequency. The second integral is divided into two parts: 

p. (cp) cos MpdCp = f^ 1 F-,(9)cos kCpdCp + f°° F x (cp)cos k9d9 (l6) 

J o o cp x 

where 9 X is chosen such that F x (<Pi) is close to zero. The first 
integral in equation (l6), being bounded, then causes no difficulty. 
Expanding the trigonometric term and retaining only the first term in 
the expansion, there results F x (cp)dcp, which is the characteristic 

area out to the point 9i* Now for large values of 9 , F;l( 9) is approx- 
imated in reference 4 by 

F i< ?) ■ “{rr? + rsip} 

where the values of p, a, and b are dependent on Mach number, and are 
given for Mq = 0, 0-5, and 0.8 in reference 4. Inserting this quantity 
in the second integral in equation (l6), we have 


f F, (9)cos k9d9 = p / 

\ 1 ff 


00 cos k9 


'9 X a 


+ 9 


. /•“ cos k9 

dcp + ** V, kf 


d9 


(17) 


Performing the indicated integrations in equation (17) , there is obtained 


-jk [a+9i ]j- 


a term, pb/a+9 x , from the second integral, and a term, -pCi- 
from the first integral. For small values of the argument, the cosine 
integral is approximated by (see ref. 12), 


Ci -j^k[a+9i ij” ~ ln-|7k[a+9 1 ]j^ 
k[a+9 x 3 — > 0 


where y is Euler's constant, 1. 78107 • 
Then, through the first order in frequency. 


Cl 


= C Lq (») - J 91 F 1 (9)*P + G(9 x ,k) 


qc/2V 0 


( 18 a) 
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where 




p. Zn-j^k [a-HPi ] j* - 


Ub 

a+<P x 


The out-of -phase moment for pitching about the leading edge follows from 
the above development, with the added result of reference it- that for 
large values of cp> 


F 3 (<p) 2 


Fi(<P) 

4 


Then, through the first order, 

- %(”) - 4® 1 FsMW - | 8(<f i ,k) ( 18 b) 


Thus, after fixing <T>i, choosing a (small) value of k, and computing 

G(<P x ,k), the finite areas corresponding to the terms - / < ^ 1 F^ (9) d<P 
m, , o 

and -f F 3 ( ( P)d < P can- be assessed in the same manner as will be done 

for the supersonic case. The advantages of such a procedure will be 
evident later. 


Correspondence between indiclal lift and moment analysis and other 
methods .- Before proceeding further with applications of the indicia! 
response method, it is appropriate to discuss the relationship of this 
approach to other widely used methods. 

Following the fundamental papers of Bryan and Routh, which intro- 
duced the basic differential equations of motion of rigid bodies and 
their stability criteria, the historical development of the theory of 
longitudinal motions of an aircraft evolved separately in two fields of 
research: dynamic stability and flutter. Workers in the dynamic stabil- 
ity field soon found that the longitudinal oscillations of a rigid air- 
craft in flight were generally of small reduced frequency. On this 
basis, the constants due to the aerodynamic properties of the airframe 
which appear in the differential equations of motion were considered to 
be independent of frequency. As a first approach to the problem of 
obtaining the necessary aerodynamic coefficients analytically, the 
instantaneous normal velocity distribution at the surface of the airfoil 
was assumed to be constant with time. The aerodynamic forces and moments 
arising from the fixed boundary conditions were then calculated using 
steady- flow theory. Later, this assumption was realized to be an over- 
simplification for the case of wing-tail combinations and an additional 
term correct to the first order in frequency was added which accounted 
for the lag in the tail pitching moment caused by the time required for 
the vorticity discharge from the wing to reach the tail (see ref. 13). 
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Since at low speeds the pitching moment of the tail far outweighs all 
other contributions, the results from steady-flow theory together with 
the term accounting for the v'orticity lag satisfactorily predicted the 
dynamic longitudinal motions of wing-tail combinations, and it waB con- 
cluded that the major aerodynamic effects had been accounted for. In 
recent years, however, numerous authors (in particular. Miles; see, e.g., 
ref. 11) have pointed out that the above-mentioned theory overlooks 
important contributions to the aerodynamic forces and moments which, 
though still within the first order in frequency approximation, arise 
from time-dependent boundary conditions and must be calculated from 
unsteady- flow theory. It has been shown by these authors that with 
proper inclusion in the equations of motion of these coefficients, the 
deterioration of damping in the short-period mode actually occurring for 
aircraft flying at speeds near the speed of sound can be successfully 
predicted. The consequences of the assumptions involved in the classical 
dynamic stability theory will be more evident from a brief review of the 
equation of motion and boundary conditions for the single-degree-of- 
freedom rotary oscillations of a rigid wing flying at constant supersonic 
speed. At the very outset, the assumption is generally made that the 
aerodynamic reactions to the motion of the airframe depend only on the 
angular position and angular velocity and not upon angular accelerations 
or higher time derivatives. The equation of motion for the change in 
pitching moment following a displacement from an equilibrium position is 
then written in the form of a power series: 


1 *rr - 

dCm"! 

^ , f Sc m 1 ac + f dCm 1 ac + 

d 2 C m 

qoSc 

_ da, J 

Id [ac/2V 0 3 i 2V 0 Id [qc/2V 0 ]J 2V 0 

_ d 2 a - 


f 1 [ac/2V 0 ] 2 t f d g Cm 1 [<x5/2Vq] 2 ; 

ls 2 [dc/2V 0 ]J 2*. I a 2 [qc/2V 0 ] -> 21 


+ 


(19) 


It should be remembered that for the rotary-oscillation case, the air- 
foil is subjected to changes in both angle of attack a and angular 
velocity q, and that these motions produce normal velocity patterns at 
the airfoil surface which are different in character. Thus, although 
for the single-degree-of-freedom case, a and q are equal, nevertheless 
their separate effects must be considered and it is therefore 

5 Cm dC™ . 

necessary to include both — r _ . =• and — — , in equation (19). 

3[qc/2V 0 ] 5 [dc/2V 0 J 

Next, if it is assumed that the moments are linearly dependent on their 
respective variables, the higher-order terms In equation (19) may be dis- 
carded and the remaining partial derivatives considered as constants for 
the given wing. There remains, therefore, a linear second-order system 
with constant coefficients. In order to calculate the coefficients 



20 


MCA EM A52L04a 


(termed stability derivatives) theoretically it became necessary, for 
lack of more refined theoretical methods, to assume that the instanta- 
neous normal velocity of the flow at the surface of the wing was fixed 

with respect to time. Thus, the partial derivative could 

S[qc/2V 0 ] 

be calculated as the pitching moment due to a, constant pitching rate, 
that is, C m{ j( cr a ) , while the derivative SCm/Sa becomes the pitching 
moment due to a constant angle of attack, that is, C^Ca). As a con- 
sequence of fixing the normal velocity pattern in time, however, it was 

Un- 
necessary to assume that the derivative r-r — . was zero. There 

& [dc/2V 0 3 

was therefore no possibility for this theory to predict the occurrence 
of dynamic instability for a wing alone, since the only damping term 
remaining is C m ^(a a )j which is always stabilizing. When the restriction 
of constant normal velocity with time is lifted, however, the assumption 
is then made that the stability derivatives in equation ( 19 ) may be calcu- 
lated separately by fixing each of the independent variables a, a, 
and q in turn with respect to time. The derivatives Cnu(a a ) 

an< * CnwCtfa) thus remain unchanged, but the derivative — ( or r. ) 

d[dc/2V 0 ] ^ 

can now be included and calculated as the pitching moment due to a con- 
stant vertical acceleration, d,V 0 . It should be emphasized that while 
CmqCtfa) an ^ £]%,(<%) calculated from Bteady-flow theory by virtue 

of the assumed Invariance with time of the normal velocity pattern, 
must be calculated from unsteady- flow theory since for constant & 
the angle of attack varies linearly with time, as does the normal veloc- 
ity of the flow at the surface. 

It is clear that since the stability derivatives in equation ( 19 ) are 
assumed to be Independent of the frequency, the result for the aerody- 
namic pitching-moment coefficient is thereby limited to one that is 
correct only to the first order in frequency. For the single-degree-of- 
freedom case, then, the in-phaBe and out-of -phase components of the 
total aerodynamic pitching -moment coefficient correct to the first order 
in frequency, become 


% ■ W»a> 

“=72 v o = c mq( 0 a) + <Wc a ) 


( 20 ) 


By comparison with the first order in frequency result from the indicial 
response analysis (eq. (15)), it is evident that the two results are 
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Identical if the quantity -/ a F 3 (<p)drp can he shown to he equivalent 

o 

to C m i(a a ) . To show this equivalence, consider a wing, initially in 
level steady flight, which is suddenly forced down with constant vertical 
acceleration &V 0 . As seen in sketch (l), the angle-of-attack variation 
in this case is a = at, where & is a constant. Then applying the 
counterpart of equation (3) for the 
pitching moment 

Cm(t) = Cm a (r)a(t-T)dT * f 


Inserting a(t-T) = a[t-T], and per- 
forming the indicated differentiation, 

(^(t) = a J C ma ( t) dT 
o 

Now replace C^Ct ) hy Cm a (t)-F 3 ( T ) o 
and let t he greater than t a . 

Then 



KO) l(t) 


Sketch (l) 


Cm(t) = a(t)Cnu(t a ) - a / a F 3 (T)dr 

o 

and nondimensionalizing, hy replacing t and t a hy ccp/2 V Q and co a /2V 0 , 
we have 

Cm(cp) = a(cp)Cm a ( 0 a ) - — / 0a F 3 (cp)dcp; q> >o a (21) 

2V 0 ° 


Thus, the pitching moment proportional to the constant vertical acceler- 
ation parameter dc/2V 0 , which is synonomous with the definition of the 
stahility derivative C^, is found to he equivalent to the pitching- 
moment contribution due to a for the first order in frequency-rotary- 
oscillation case.® Therefore, the results of the indicial response 
method, when reduced to the first order in frequency for supersonic 
speeds, are identical to the results from the familiar first-order theory 
used in dynamic stahility work. 4 

®By the same procedure, the stahility derivative Ci^Oe,) can he shown 
to he equivalent to the term - f^ a F 1 (cp)dcp. 

4 Notice in equation (21) that if a a is replaced hy infinity, the results 
apply to subsonic speeds. For the two-dimensional wing, the analogy 

between and - f F 3 (cp) dq> then gives only the previously mentioned 

singularity at infinity as k — > 0. If the area corresponding to 
J o F 3 (cp)d9 were finite, however, the analogy would he equally useful 
for subsonic as well as supersonic speeds. 
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Workers in the field of flutter, who were concerned with frequencies 
many times those encountered in dynamic stability analyses, required 
theoretical information showing the behavior of the forces and moments 
as affected by the frequency of oscillation, and therefore discarded the 
first-order theory for more precise methods of analysis. One of the 
most useful of these has been the "oscillating potential" theory, which 
is based on periodic solutions of the time -dependent linearized equation 
of compressible potential flow. The in-pha8e and out-of -phase lift and 
moments are thereby determined, generally as functions of powers of the 
reduced frequency, aspect ratio, Mach number, and position of the axis 
of rotation. The application of this method, which developed primarily 
as a result of Theodorsen's work for incompressible flow (ref. l4) lias 
recently produced a number of useful papers covering a wide variety of 
wings at supersonic speeds (see, e.g., refs. 15 to 19) • It has been 
shown by a number of authors, in particular, Garrick, in reference 7> 
that through the use of superposition methods the results for the aero- 
dynamic coefficients obtained from the oscillating potential theory are 
wholly compatible with those of the indicial response method for har- 
monic motions (eqs. (12) and (l4)). 

Thus the indicial response method embraces both the first-order 
theory of dynamic stability and the oscillating potential theory, and, 
in effect, bridges the gap between the fields "of dynamic stability and 
flutter. Furthermore, the indicial response method overcomes the dis- 
advantages of the two methods described above since, unlike the first- 
order theory, the effects of frequency on the aerodynamic coefficients 
can be determined, and, unlike the oscillating potential theory, the 
method can be applied easily to the study of arbitrary motions. Finally, 
the indicial response method represents a fundamental approach to the 
problem of unsteady flow, and affords valuable insight into the physical 
nature of the aerodynamic phenomena taking place. 

Physical concepts relating to the indicial loading .- It has been 
shown that for even small frequencies, the pitching moment of an airfoil 
in harmonic rotary motion can lag behind the angle of attack of the air- 
foil. The magnitude of the lag depends on the character of the indicial 
response to a step change in angle of attack. It is therefore of 
interest at this time to re-examine the physical nature of the flow 
that contributes this lag. 

Consider first the lift and moment at the instant the angle of 
attack changes, assuming that previous to time zero the wing has been 
flying a level path at zero angle of attack. At t = 0, the wing begins 
to sink, without pitching, With constant downward velocity aV 0 while 
maintaining its forward velocity. The angle of attack therefore changes 
discontinuously from zero to a constant a. At the same instant, the 
step change in angle of attack causes the emission of a compression wave 
from each point on the lower surface of the wing and expansion waves 



NACA RM A52L04a 


23 


from points on the upper surface. In the infinitesimal time during which 
the starting action occurs, each section of the wing experiences the 
same impulsive force, and by equating the impulse to the momentum trans- 
mitted to the mass of fluid affected by the starting waves, the starting 
lift coefficient can easily be derived as iki/Mo (see ref. 3)* During 
the infinitesimal starting time, the pressure disturbances from the edges 
of the airfoil, propagated at the speed of sound, travel an insignificant 
distance and do not influence the remainder of the airfoil. The lift 
coefficient is therefore independent of the wing plan form. This remark- 
ably simple result for the starting lift coefficient, which is valid for 
both subsonic and supersonic speeds, is thus dependent solely on the 
flight Mach number. The starting pitching moment follows directly from 
the above result, since by virtue of the uniformity of loading the aero- 
dynamic center is located at the wing centroid of area. 

For values of time greater than zero, however, the situation differs 
radically for the supersonic and subsonic speed ranges. Consider first 
the supersonic case. As time passes, the spherical sound waves emitted 
at t = 0 grow in size with radius ad . The wing, however, is moving 
forward at a faster rate than the rate of growth of the starting sound 
waves and thus begins to emerge from the influence of these waves. This 



Sketch (m) 


At t = 0, the starting waves just cover the wing and the loading is 
uniform as described previously. At t = t x , the starting waves have 
grown in radius and the wing has begun to emerge from their influence. 

On that portion of the wing which has emerged, region (l) in sketch (m) , 
the loading has already reached its steady-state value. Notice that in 
this region the characteristic tip Mach cone has already formed. On the 
portion of the wing uninfluenced by the starting waves from the edges, 
region (3) in sketch (m), the loading is still uniform as at t = 0. 

In regions (2) and ( 4) , the loading is Influenced by the starting waves 
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from the leading and side edges, and in these regions is thus different 
from the loading in either regions (l) or (3)* As time increases still 
further, the uniform starting load quickly disappears as the sound waves 
from the leading edge grow in size and as the wing moves forward. 

Finally, at time ta, the envelope of the starting waves from the lead- 
ing edge is coincident with the trailing edge of the wing, and the 

steady-state loading corresponding to the new angle of attack a has 

been completely established over the wing. 

The above relationships can be shown more clearly for the entire 
time interval zero to t a for a two-dimensional wing by plotting as a 
function of time the position of the wing leading and trailing edges 
and position of the envelopes of the sound waves which emanate from the 
leading and trailing edges at t = 0. Sucha plot is shown in sketch (n). 



It is clear that at t = ti the regions of the wing (l), (2), 
and (3) correspond to the same regions at t = tj. for the wing shown 
in sketch (m) . For t = 0 and in region (3) the loading is uniform and 
is given by ka/Mo. For t>t a and in region (l) the wing has out- 
stripped the starting waves from the leading edge and has attained its 
steady-state loading. For t<t a the chordwise loading is composed of 
combinations of the loading in each of the three regions as shown in 
sketch (o) . 5 


5 The reader will note the similarity between sketch (n) and sketches 
depicting the boundary conditions for three-dimensional wings in 
steady supersonic flow. Many researchers have pointed out the analogy 
and it has been used to calculate the pressure over a wing impulsively 
starting from rest (refs. 3 and 4). 
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It may be easily verified that t a , the time required for the wing to 
attain its steady-state loading, is given by the following relationship: 


ta — 


c qMq 

V Q [Mq-csc v] 


<£ 


( 22 ) 


Notice that for wings having straight or swept-forward trailing edges 
and straight or swept-back leading edges, equation (22) reduces to 


c qMq 

V 0 [Mo-1 } 


(23) 


In terms of the number of half M.A.C. lengths of travel, it thus appears 
that for wings having straight or swept-forward trailing edges and 
straight or swept-back leading edges, 


- tzO® t - 2c O Mp 
C a " c [Mo-1] 


(24) 


whereas for wings having swept-back trailing edges and straight or 
swept-back leading edges. 


°a = 


2c 0 Mq ^ 2cq 

7C* Mo sicsc V + cos V 

c [Mq-csc VJ ^ b 


cfa — 2 Mq 


Mq 




4c c 


[Mo 2 -i3 


1 

Mq^csc V + — - 


1(25) 


COS V 


Mq 2 - 1 


The second of equations (25) applies to that range of Mach numbers for 
which the trace of the starting sound wave from the apex is not tangent 
to the trailing edge at t = t a . 


Now consider the subsonic case. Here the situation Is more compli- 
cated In that since 'the starting sound waves travel faster than the 
wing, the wing never escapes their influence. Furthermore , the vortic- 
ity shed by the wing at t = 0 can also influence the loading on the 
wing since the disturbances created by the shed vorticity travel forward 
at a faster rate than the wing. For these reasons, the indicial loading 
at subsonic speeds approaches Its steady-state distribution asymptoti- 
cally with time. The situation for subsonic speeds will be more clearly 
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and in the region corresponding to region (3) of sketch (q), the loading 
is uniform and equal to 4a,/Mo, as in the supersonic case. 


Damping in Pitch of Low-Aspect-Ratio Wings 


Previously (eq. (20)), it was shown that for single-degree-of- 
freedom, slow- frequency, pitching oscillations of a wing, the principal 
parameter contributing to the damping of the motion is the damping coef- 
ficient + Cmrj - This result, however, is not directly applicable to 

the case of an aircraft in flight, since generally additional damping is 
provided by virtue of the fact that the aircraft experiences harmonic 
vertical translatory oscillations as well as the rotary oscillations. 

For the s low-frequency case, however, the effect of the second degree of 
freedom on the damping can be determined by use of the relationships 
given in reference 10. It is further shown in reference 10 that although 
the effect of the translatory oscillation is usually to increase the 
total damping, nevertheless, the parameter of primary importance in 
determining the magnitude and duration of the motion remains the damping- 
in-pitch coefficient 0^ + Cm- . The remainder of this section is 
therefore concerned witha study of the effect on this parameter of 
certain important variables. In particular, the effect of the position 
of the center of gravity, and the effects of Mach number, aspect ratio, 
plan-form shape, frequency, and thickness will be examined, principally 
by inspection of the indicial lift and moment responses to a change in 
angle of attack. 

Effect of static margin .- From the previous discussion, it will be 
remembered that at supersonic speeds the stability derivatives cr a ) 

and C™. ( ca) were shown to be equivalent to the indicial lift and moment 

expressions, - f F 1 (cp)dxp and - / F s (cp)dcp. Furthermore, it was shown 

that these were the contributions which could cause the total lift and 
moment during the short-period oscillation to lag behind the angle of 
attack. Hence, by inspection of equation (15), it is evident that since 
Cmq( cr a ) is always stabilizing, when C ma (a a )is negative (corresponding 
to a statically stable condition) the possibility of dynamic instability 
in the form of negatively damped rotary oscillations arises when 
Cjn^Oa), the shaded area in sketch (k-ii), is larger than C 1D g(o’a)* 

Now since the normal velocity at the surface of the wing due no the 
instantaneous angle of attack is constant over the wing, the lift 
derivative Cl^ which arises from this boundary condition is independ- 
ent of the axis position and Cm^ will therefore vary linearly with 
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axis position. This variation is illustrated in sketch (s), -where Axo 
represents the distance of the center of gravity from the aerodynamic 
center, measured positive forward of the aerodynamic center. 



The parameter Ciq, on the other hand, is a direct function of the 
axis position, since it arises from a normal velocity distribution that 
varies directly as the distance from the axis. The moment coeffi- 
cient will therefore vary as the second power of the axis posi- 

tion, and describes the parabolic shape shown alBO in sketch (s). It 
is evident from sketch (s) that the sum of Cjq^ and will be a 

minimum at some value of the static margin, and that the sign 
of Cn^ + Cmd, at that point determines whether or not a region of axis 
positions will exist over which the wing can experience negatively 
damped rotary oscillations. These qualitative statements may be written 
explicitly by considering the equation for the damping in pitch about an 
arbitrary axis. 


c mq + Cmi 



(26) 


where again AXq refers to the distance of the center of gravity frc-i 
the aerodynamic center, and the subscripted terms are referred to an 
axis through the aerodynamic center. 
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Taking the derivative of equation (26) with respect to Axq/c, and 
setting the result equal to zero, there is obtained the axis location 
at which, the damping in pitch is a minimum 


d [Cmq+Cm^ ] 
d[AXo/c] 



I 4*0 

k -T CL a = 0 




(27) 


When equation (27) is inserted into equation (26), the minimum value of 
the damping in pitch is given as 


Cmq+Cmcx 


min. 





and hence, a region of instability will exist if 


( 28 ) 


{ 




(29) 


If equation (29) is greater than zero, the boundaries of the region of 
axis positions over which instability is possible are of course given by 
setting equation (26) equal to zero and solving for Axq/c. 


AXq 

c 





Notice In both sketch (s) and equation (30) that for a given Mach number 
there will be two axis positions at which the damping in pitch vanishes. 
Then if the above procedure is carried out for a series of Mach numbers, 
one may trace out a curve as shown in sketch (t) which forms the locus 
of Mach numbers and axis positions at which the damping in pitch is zero. 
This locus thus delineates the regions of Mach number and axis position 
for which dynamic Instability iB and is not theoretically possible. 

Such loci, covering a wide variety of wings at supersonic speeds, have 
been presented by a number of authors. Watkins, for example, presents 
supersonic boundary curves for rectangular and triangular wings in refer- 
ences 16 and 17. At subsonic speeds. Miles' reduction of Possio’s 
development to first order in frequency (ref. 11) can be used to form a 
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stability boundary curve for the two- 
dimensional wing for a given (small) 
reduced frequency. Such a curve is pre- 
sented for the entire Mach number range 
in figure 1, where, here + Xq is the 
distance of the axis of rotation behind 
the leading edge and k, the reduced fre- 
quency, is 0.011 for subsonic speeds and 
approaches zero for supersonic speeds. 

Notice in figure 1 that at both subsonic 
and supersonic speeds, the range of Mach 
numbers over which dynamic instability is 
possible is largest for center-of-gravity 
positions forward of the aerodynamic 
center. Further, the largest range of 
axis positions over which dynamic insta- 
bility is possible occurs near Mq = 1. 

Both of these characteristics are true as well for three 
wings at supersonic speeds (see, e.g., refs. 16 to 19 ). 



•dimensional 


Effect of Mach number .- Next consider the effect of Mach number on 
the damping in pitch of a two-dimensional wing with axis at the leading 
edge. The variation with Mach number of the indicial pitching-moment 
response to a change in a will first be examined, using the information 
given in the previous sections and the indicial curves given in refer- 
ence A. At supersonic speeds, the manner in which - / o aa F 3 (cp)dCp, the 


area corresponding to Is affected can then be assessed and com- 

pared with Cmq ( era) • At subsonic speeds, UBe is made of equation (l8b) . 
It is evident in equation (l8b) that by fixing k and choosing <Pi, 

such that the quantity — G(cp 1 ,k) is the same at each Mach number, one 

k cp. 

is free to compare finite areas - J F 3 (cp)dcp on an equivalent basiB. 


As has been mentioned previously, the starting lift, at any Mach 
number, is ka/Mo and is concentrated at the midchord. At Mq = 0, 
therefore, there is an initial infinite pulse in the pitching moment 
about an axis coincident with the leading edge after which the indicial 
curve drops to jt/k and begins to grow asympotically toward its steady- 
state value jt/2. At low subsonic Mach numbers other than zero, the 
initial pitching moment is less than infinite but very large, being 
2/Mo, and then falls before growing toward the steady value jt/20. As 
the Mach number increases toward 1.0, the starting pitching moment 
falls while the asymptotic value grows, until at Mq = 1*0 the indicial 
curve becomes unbounded in asymptotic moment. 
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Ab seen in sketch (u), the effect of increasing the Mach number at 
subsonic speeds is therefore to increase rapidly the area corresponding 

to the destabilizing moment contribution, - F_(cp)d<p. 



Sketch (u) 

In sketch (u), k was chosen to be 0.011 and" the values of were 

picked such that - i G(q> 1 ,k) was +4.88 at each Mach number. In the 
following discussion, the damping moment - F 3 (cp)dtp + 4.88 will be 
referred to as Cm^ for convenience. ° 

At supersonic speeds, the initial value of the pitching moment 2 /Mq 
continues to drop with increasing Mach number, but here the steady-state 
pitching moment also begins to fall and at a faster rate than the start- 
ing moment, being 2/0. Even more important, as the Mach number increases, 
the number of half-chord lengths traveled to: reach steady state decreases 
rapidly, being 22, for example, at Mo = 1.1, as compared to 4 at Mq = 2. 
As seen in sketch (v), the area representing therefore shrinks 

rapidly with increasing supersonic Mach number and becomes relatively 
unimportant at Mach numbers greater than 2. The trend of with 

Mach number through the range 0 <Mq< 2 is more clearly evident in 
sketch (w) . It is seen that is positive, or destabilizing, 

throughout the Mach number range and that its effect is most important 
at Mach numbers near 1.0. AIbo shown plotted for comparison in sketch (w) 
is the variation of C m q with Mach number. When the parameters are 
added, it is evident that the damping moment for the two- 

dimensional wing with axis at the leading edge is destabilizing in the 
Mach number range 0 < Mq < 1 . 4l4 . 
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Sketch, (v) 


Effect of aspect 
ratio . - To illustrate the 
effect of aspect ratio, 
it is convenient to 
compare the supersonic 
damping-in-pitch char- 
acteristics of a group of 
triangular wings having 
subsonic leading edges. 
The wings are of equal 
area and differ only in 
aspect ratio. As was 
done previously, the 
indie ial lift responses 
to a change in angle of 
attack will first he 
examined. The effect of 
aspect ratio on the char- 
acteristic area repre- 
senting can then he 

assessed. 



Sketch (w) 


As has heen mentioned previously, the starting lift coefficient 
after a step change in a is independent of aspect ratio and is there- 
fore equal to 4a,/Mo for each wing. The parameter a a , the number of 
half M.A.C. lengths required to reach steady state, is also the same for 
each wing, being a function only of Mach number. The steady-state lift 
coefficient, on the other hand, is a function of aspect ratio and 
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±n C «!2+ e £ f s x t5ae as P ect ratio ia reduced, 
in sketch (x), as the aspect ratio becomes 


Thus, as shown schematically 
smaller, the characteristic 



Sketch (x) 

isltsts ^ b d 8 ore : t T rapi<iiy - 0 ^ -u«t 

state lift coefficient ? P° 8lt ^e since the area below the steady- 
It is evident, thllSlrl St Til aSf'f*' 1 f ° r iy the above, 

stabilizing effect on the damping lIliS LIS*? ^ “ hlBhljr 

Of the static margin the develnJLnt f for P ositive values 

is possible only when cv -ton , , f a d - es ' fca hilizing damping moment 
8 ketch (y) S for ' ThU re8ult is in 

the dam JE, If "* *W- 2 - 

Since for triangular wings the lift due to Til ollllZZTT 1 , 
^ equal to - | C^. The variation of with aspect rati! 

in Ketch’wf wlo'pStSfrsSlIfcyJ^rtS^SiS sh ° m 


sketch Fv^ !' DBOnlc doling edges. The curves drawn in 

tTSicSe^I. true eh ^ s ’ nre intended oSy 

The ISIt lallltS If r ^acteristio area with aspect ratil. 
the results of reference M. “ “ aspeot rati0 oan he computed from 
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with aspect ratio (ref. 20). It is apparent that although Cm^ becomes 
more stabilizing with increasing aspect ratio, the destabilizing effect 
of Cm^ predominates, and the trend of the net damping moment is seen 
to become highly destabilizing as the aspect ratio is increased. 


By the same reasoning, "the variation with aspect ratio of the 
damping moments of other types of wings can be shown to be similar (see, 
e.g., refs. 11, 21, 22, and 23). A notable exception, however, is the 
triangular wing with supersonic leading edges, whose damping in pitch 
has been shown to be independent of aspect ratio (see refs. l8 and 19 )* 
This characteristic may be anticipated from a study of the indicial 
response curves, since not only are the initial pitching moment Cjn^o) 
and the half M.A.C. lengths traveled to reach steady state (cr a ) inde- 
pendent of aspect ratio, but, unlike the sub sonic -edged triangular wing, 
the steady-state pitching moment Cn^cTa) is also independent of aspect 

ratio, being - — — — . Inspection of the results of reference 4 then 

p L c J 

reveals that the indicial variation Cm^cp) between zero and o a and 
the steady-state parameter C n!( ^(cr a ) are likewise independent of aspect 
ratio . " 


Effect of plan- form Bhape .- Next, consider the effect of plan-form 
shape on the supersonic damping in pitch of a group of wings having the 
same aspect ratio. In order to make use of readily available theoretical 
results and yet compare these results with the experimental data to be 
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given later, three wings of aspect ratio 3 are chosen, having the trian- 
gular, swept, and rectangular plan forms shown in sketch (z). 



Sketch (z) 

As has been shown in the section entitled "Effect of static margin, 11 
the tendencies of the wings toward dynamic instability can be compared 
comprehensively by plotting their stability boundaries. For this com- 
parison, then, use is made of equation (30). The stability derivatives 
which appear in equation (30) were computed from the theoretical results 
of references 11, 20, 22, and 23- Results of. these calculations are 
shown in figure 2, wherein the stability boundaries for the three wings 
are shown as a function of axis position and Mach number. (Note that 
the axis position for each wing is measured as the distance from the 
leading edge of the M.A.C. of the wing, and that the dimensions are non- 
dimensionalized on an equivalent basis by referring them to the M.A.C. 
of the triangular wing.). It is clear from inspection of figure 2 that 
at any Mach number the triangular wing has the smallest range of axis 
positions over which dynamic instability is possible and the rectangular 
wing, the largest. 

The differences in the damping characteristics of the triangular 
and rectangular wings will be more clearly understood by a qualitative 
study of their indicial lift responses for a Mach number of 1.2, and an 
examination of the distribution of loading due to d for the two wings. 
Consider first the indicial lift responses. 

Again, the starting lift coefficient is independent of plan -form 
shape and is ta/Mo for each wing. For the rectangular wing, the 
lift drops abruptly after time. zero due to the loss in lift in the 
regions of the wing influenced by the formation, of the tip Mach cones 
and the starting waves from the side edges (see sketch (m)). Then as 
the wing begins to emerge from the influence of the starting waves, the 
lift begins to recover, rises toward its steady-state value (given 
by eq. 6.3-2 of ref. 11), and attains this value after 12 half-chord 
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lengths of travel (eq. (24)). The variation is shown in reference 4 and 
is redrawn in sketch (aa). As mentioned previously, theoretical indicial 



Sketch (aa) 

lift results have not yet been developed for triangular wings having 
subsonic leading edges. However, the variation shown in sketch (aa) is 
considered to be a reasonable estimate of the true shape, being based on 
knowledge of the steady- state lift (ref. 20 ), the time to reach steady 
state (eq. (24) ) , and the assumption that the shape of the variation 
would be similar to that of the wide triangular wing (ref. 4), The 
curve was adjusted within the known end points until the area correspond- 
ing to Cj^ agreed with that given for this parameter in reference 20. 

It is evident from examination of sketch (aa) that because of the initial 
loss in lift and the larger steady-state lift for the rectangular wing, 
Ci^ for this wing is significantly more negative than that for the tri- 
angular wing. Next, it is shown in references 11 and 21 that with the 
exception of regions Influenced by tip Mach cones, the loading due to a 
(for cp > a a ) for wings having swept-back leading edges increases linearly 
from zero along rays from the apexj whereas for rectangular wings the 
loading due to d Increases linearly from zero along chord lines. These 


characteristics place the center of loading due to 

2 


Ot* 0/fc * Cq 


for the 


triangular wing and approximately — c Q for the rectangular wing. 7 

Then for an axis of rotation passing through the aerodynamic centers of 
the wings, the moment arm for the lift due to d for the triangular 


7 Due to the influence of the tip Mach cones, the center of loading due 
to d is shifted forward somewhat from the position it has for the 
two-dimensional wing. Calculations for the A = 3 rectangular wing 
at Mo = 1.2 show that the center of loading is at 0.605 c Q . 
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wing is 



whereas for the rectangular wing it is approxi- 


mately 



Thus, not only is the negative out-of -phase lift contri- 


bution Cl^ for the rectangular wing significantly larger than that for 

the triangular wing, but the destabilizing damping moment - $ Cjr is 

c 

larger yet, due to the larger moment arm. Calculations for the steady 
pitching parameter for an axis through the aerodynamic center 

(refs. 11 and 20 ) then reveal that for the triangular wing is more 

negative than for the rectangular wing. The net result is therefore a 
considerably larger damping moment for the triangular wing than for the 
rectangular wing. The result of this comparison, however, should not be 
interpreted as a recommendation that the triangular rather than the rec- 
tangular wing be used on aircraft from a dynamic stability standpoint. 

To obtain adequate static stability, the rectangular wing would generally 
be employed in combination with a tail surface, whereas the triangular 
plan form may be sufficiently airworthy without the use of a tail. The 
addition of a tail surface in effect reduces the aspect ratio of the 
rectangular wing, which reduction, as noted previously, haB a highly 
stabilizing effect on the damping in pitch. The tailless triangular wing 
may therefore experience more difficulty at Mach numbers near 1.0 than a 
rectangular wing-tail combination. 


Effect of frequency .- The previous discussion has been restricted 
to the analysis- of a harmonic motion that is of vanishingly small fre- 
quency. This limitation. arose as a consequence of discarding all but 
first order in frequency terms in the expansions of equations ( 12 ) 
and (l4). The question arises: when the frequency can no longer be con- 
sidered small, what effect has the frequency on the damping in pitch? 


Previously, the trigonometric terms in equation (l4) were expanded 
and, assuming k to be very small, terms of order k 2 and higher were 
eliminated. It was then found that the loss in damping from that pro- 
vided by the steady damping parameter Cmq^Ca) was associated with the 
destabilizing moment contribution corresponding to the term 

- f Q a F 3 ( <P) d<P. Now, however, we discard the restriction of small k 
and perform graphically the integrations evident in equation ( l4) for 
supersonic speeds for several values of k. The procedure is Indicated 
in sketch (bb). It is apparent from sketch (bb) that the effect of 
increasing k is to reduce the area corresponding to the destabilizing 

moment contribution - J a E Q ( 9 ) cos kcpdcp. There appears another 


8 Again, due to the influence of tip Mach cones, the aerodynamic center is 
shifted forward from 75 c Q to 0.W+3 c Q . The moment arm is therefore 

O.163 c Q . . 
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Sketch, (bb) 

destabilizing contribution, -k P 4 (q>) sin k<pdq>, but quite evidently 
its effect is sma ll compared to the reduction in the term 
. / aa F s (9) cos kcpdcp. Notice further in sketch (bb) that the effect of 
increasing k becomes of marked importance -when the half -period of 
oscillation is the same order of magnitude as the time for the indicial 
response to reach steady state. As shown in sketch (bb) for the fre- 
quency k 2 , thq area^b^then begins to subtract from(a^, so that the 
destabilizing contribution - F g (<P) cos ktpdcp can be very small. We 
may therefore expect that increasing the frequency of oscillation has a 
stabilizing effect on the damping in pitch. This conjecture is sub- 
stantiated in figures 3 and b, where the supersonic stability boundary 
curves for aspect ratio k triangular and rectangular wings are shown 
plotted for various reduced frequencies. These curves were obtained 
from calculations based on the results of references 16 and IT* Notice 
that for both wings the region of possible instability is diminished as 
the frequency is increased. 

From the results of the analysis for supersonic speeds, we may 
further expect that the stabilizing effect of increasing the frequency 
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will be of even more importance at subsonic speeds, for here the indieial 
variation F s (<3>) dies out at infinity. The half -period of oscillation 
is therefore always smaller than the time to reach steady state. 

QQ 

The situation for the term - f F a ( c P) cos kcpd<p is shown in sketch (cc). 



Sketch (cc) 

It is evident in sketch (cc) that the destabilizing moment 

- J F a (q>) cos k<Pd<P diminishes rapidly as the frequency is increased. 
The effect of this reduction on the damping in pitch can be illustrated 
by plotting the subsonic damping-moment coefficient against reduced 
frequency for the two-dimensional wing (with axis at the leading edge) 
for Mach numbers 0, 0.5> and 0.8. The results, which were obtained 
from reference k, are shown in figure 5 . It is seen in figure 5 that 
the large destabilizing effects of the moment contributions 

- J P 8 (<P) cos kcpd9 and -k f o F 4 ( C P) sin k<Pd9 are confined to a rela- 
tively narrow range of reduced frequencies. Notice further in figure 5 
that the range of frequencies for which instability is possible is small 
at Mq = 0 (0<k<0.o4) and grows with increasing Mach number. This is 
believed to be the primary reason why unsteady lift effects were found 
to be unimportant at low speeds but are of great importance at speeds 
near the sonic speed. 
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Effect of thickness .- In this section, the concepts regarding the 
indicial functions discussed in previous paragraphs will he used to 
examine the effect of profile thickness on the damping in pitch of a 
two-dimensional wing flying at a supersonic Mach number - 

It will be recalled that for low frequencies the damp ing in pitch 
may be considered to be the sum of the damping moment due to steady 
pitching, Cm^cTa}, and the area representing on the plot of the 

pitching-moment response to a step change in angle of attack 
(sketch (k-ii)). The airfoil thickness, of course, influences both 
these parameters. For the purpose of this analysis, it will be con- 
sidered sufficiently accurate to study the effects of thi ckn ess only on 
the steady parameter C^Ca) and on the end values C^Co), C^Oa) 
and 0 a of the indicial response curve. It will then be assumed that 
shape of the indicial response curve is not appreciably altered by thick- 
ness effects, whereupon the major effect of thickness on can be 

assessed by adjusting the response curve given by the linearized theory 
to fit the corrected values of C^o), C^Ca) and Sa- 
lt has been shown by Busemann (see, e.g., ref. 24) that the loading 
at an element of a lifting airfoil may be represented by 

= C 1 [a u -0 l +2a] + C 2 jcr u 2 ~cr i s +2a [a u + 0 Z ] j- (3l) 

where 

J m/-i 

E 7+1] Mp 4 - 4[Mo g -l] 

2[Mo 2 -1] 2 

and 0 U and o z are the local slopes of the upper and lower surfaces, 
respectively, measured with respect to the chord line. For the special 
case of an airfoil that is symmetrical about the chord line, 

0 U equals 0 Z and equation (31) reduces to the following: 

= 20^! + 4oC 2 a u (32) 

For the sake of simplicity, in what follows equation (32) will be used 
for the local loading rather than equation ( 31 ) • The results will there- 
fore be valid only for airfoils that are symmetrical about the chord line. 
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The lift on the airfoil due to angle of attack is obtained simply by 
integration of equation ( 32 ) across the chord. 


Cp/ CT a) 


-i- f C °^ dx 

c O O 'lo 


CLo (, a) = ^ 


2 Ci 



(33) 


Similarly, the moment due to angle of attack, measured about the leading 
edge, is 


<vK) 




Ap 

9o 


x dx 


Cmo,' ( cr a ) 


Cm 1 ( g a) 
a 



m 


When both the leading and trailing edges are on the chord line, equa- 
tions ( 33 ) and ( 3 k) reduce to the following: 


Ci Ja a ) = 2C X 


Cmo , 1 (°a) C a 



(35) 


where ft is the enclosed area of the airfoil section. 


The local loading given by equation (32) is still applicable for 
the steady pitching case with this added provision: that now a repre- 
sents the local angle of attack. For an airfoil pitching about its 
leading edge, the variation of a with chordwise distance is given 
by a = qx/V Q , where q is the pitching velocity. Then again, for the 
lift due to pitching. 


C L ’ (d a ) 


= i/ C ° £Ecbc 
C 0 O 3o 


<V(Os.) 


CL'(da) ^ . 8C 2 r c Q 

= + — 5 J 

qco/2Vo c o 0 


c^x dx 


(36) 
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Likewise, the pitch ing moment due to pitching about the leading edge is 


Cm' ( ff a) 


-2* ; Co & * 

Co 2 ° 


dx 


Cmq’ (cr a ) 


Cm' ( g a) 

<lc 0 /2V 0 




(37) 


The primes on the parameters indicate that the moments are measured about 
and the wing is pitching about the leading edge. 

In order to illustrate the result given by equation (37)# the 
damping moment due to steady pitching has been computed for a family of 
airfoils having symmetrical biconvex parabolic arc sections. The 
results are presented in figure 6. The dashed line in figure 6 repre- 
sents the theoretical locus of Mach number and leading-edge angle for 
bow-wave detachment (ref. 24). Since the requirements of the second- 
order theory are violated to the left of this line, the damping curves 
there are not rigorously correctj and therefore are terminated a short 
distance beyond the line. 

To determine the effect of thickness on the parameter, Cm^' , it is 
assumed that the shape of the indicial pitching-moment variation is not 
significantly affected by thickness. When this assumption is made, it 
becomes necessary only to correct the initial and final ordinates of the 
indicial curve Cm^' (o) and Cmo,' (era) and the number of half-chord lengths 
traveled to reach steady state, cr a . The response curve 1 b adjusted to 
fit these new end values, so that the area representing the cor- 
rected can then be determined. 

The effect of thickness on the final ordinate of the indicial curve, 
Cn^' (oa)# has already been determined (eq. (34)). NoV consider the 
starting lift and moment. 

Assume first a thin flat wing to be flying at zero angle of attack 
in a uniform stream of density p Q . At time zero, it Btarts to sink 
with downward velocity V G a. The impulsive start causes a plane com- 
pression wave to be emitted from the lower surface at the speed of sound, 
which is constant throughout the stream. (The pressure and therefore 
the temperature is sensibly constant.) At the same instant, a plane 
expansion wave is emitted from the upper surface and travels at the same 
velocity in the opposite direction. At the end of an infinitesimal time 
period At, the wing has sunk a distance V 0 oAt. The sound waves have 
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moved a distance aoAt in both directions. The mass of fluid affected 
is 2p 0 c 0 a 0 At. Then, from Newton’s second law, one may write_ 

(fu”f z) At = 2p 0 c 0 aoAt V Q a 

where *VL and f 2 are the forces per unit span on the upper and lower 
surfaces, respectively. Converting to coefficient form, it then appears 
that 


C L (o) _ f u -f Z = _4 
a \ PoV^o* “o 

For the wing of finite thickness for which the pressure varies along the 
chord, a different result is expected. The wing is assumed to have been 
flying at zero angle of attack for a period of time long enough for the 
flow to have established itself around the surface. In this case, the 
density, temperature, and therefore the speed of sound, vary along the 
chord due to the curvature of the surface. At time zero the wing again 
starts to sink, but the compression and expansion waves are propagated 
at different rates at various chordwlse positions so that 


At [f u -f J(x) = 2 (x) (x) Poho^t V D a 

Po a o 

and the starting lift becomes 

Cl( - ) - s — iL- f c ° JL JL dx (38) 

a ^o c o 0 Po a o 


From reference 24, the variation of p/p 0 a / a o with Cp is found to be 


_P_ ja_ 
Po a o 




(39) 


where C P is the pressure coefficient at a point on the surface when 
the airfoil is at zero angle of attack. Again, from reference 24, 


C P " 


P-Pq 

% 


^1°U + ^£Pu 




where p and p 0 are, respectively, the pressure at a point on the sur- 
face of the airfoil and the free-stream static pressure. 
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For tiie symmetrical biconvex parabolic arc airfoil considered 
earlier, the local slope o u is given by 


a u = 25 


1 - — 

=o J 


( 41 ) 


where 8 is the thickness ratio. 


Then expanding equation (39) in a binomial series, retaining terms 
through the second order in a u , and performing the integration indicated 
by equation ( 38 ), there results for the starting lift coefficient 


c La(°) 


~ -jl+0 . 04 Mq 2 5 2 [20C 2 -M 0 2 C 1 2 ] j 


(te) 


A similar procedure can be carried out for the starting pitching-moment 
coefficient, measured about the leading edge 


c ffla' (°) 


4 r c ° o a , 
J ±. x dx 

MqC 0 2 ° Po a o 


Cn^'Co) = - -~ |l+0.04Mo 2 [20G 2 6 2 -M o £ C 1 £ & 2 -10C 1 5 ] j (43) 

It is necessary also to consider the effect of thickness on the time 
required for the loading on the airfoil subjected to a step change in 
angle of attack to reach steady state. As was seen in sketches (n) 
and (p) , the time to reach steady state is exactly equivalent to the 
time required for the wing to escape the influence of the sound waves 
emitted from the leading edge at the start of the motion. For purposes 
of analysis, the situation is reversed by requiring that the wing remain 
stationary in a fluid moving with supersonic velocity V D , and causing 
a sound wave to be emitted at the leading edge of the wing. To an 
observer standing on the upper surface of the airfoil, there will appear 
to be two wave fronts whose velocities tangent to the surface are 
[V+a] and [V-a], where V and a are the local stream velocity and local 
speed of sound, respectively. Then the time required for the slower 
wave front to clear the airfoil is given by 


*a = 




dx 

[V-a] cos[tan _ 1 a u ] 


Mg p °° Jl+a u S dx 
v o J 0 [V/a 0 -a/a 0 I 


m 
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It ie therefore necessary to obtain V/a Q and a/a Q as functions of 
ehordwise location x. Again, in reference 24 these quantities are 
found to be functions of the local pressure coefficient Cp. 


7-1 1 



Expanding equations ( 45) through the second power of 
through the second power of cr u , there is obtained 


- 


J^o 

V q [Mo 


^c 0 [1 + | a u 2 J dx 

l-^Cp+^Cp 2 


Cp and 


05) 



(46) 


where 

0.50Mo + O.IOMq 2 

K = ' 1 ' ' 

Mq - 1 


K 


2 


0 . 12^Mq [Mq 2 -! ] + 0.03Mq 4 
Mq - 1 


The pressure coefficient Cp is now assumed to be that existing at 
a point on the surface when the airfoil is at zero angle of attack. This 
assumption, of course, introduces some error since the angle of attack 
has changed discontinuously from zero to a at t = 0. However, if we 
require that a be less than 5 2 , then the error introduced through 
neglect of the effect of the change in a on the pressure coefficient 
is negligible. Then, as in the starting lift problem, 


C p = C 1 c u + C e a u a (47) 

Inserting equation (47) into equation (46), dividing numerator by denom- 
inator, and again retaining terms through the second power in 0 U , there 
results 


- HZ 
VMq-13 


f 


1+C l K l ff U 



k iC 2 +E | V 2 ~k 2 3c 1 2 + | | 




dx (48) 
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For the symmetrical biconvex parabolic arc airfoil, o u is replaced by 
equation ( 41) , whereupon for this airfoil the time to reach steady 
state becomes 


t a 


— ^ [ 1 + 6 ] 

V 0 [Mo-l] 


(49) 


where 


£ = 



Alternatively, in terms of the number of half-chord lengths of travel. 


CT a 



2Mp 

Mq-1 


El+13 


(50) 


The effect of thickness on the damping parameter 0^* can now be 
estimated approxi ma tely by adjusting the indicial pitching-moment curve 
obtained from the linearized theory to the end values given by equa- 
tions (34), (43), and ( 50 ), and then measuring the area corresponding to 
the "corrected" Cm^. The process is indicated in sketch (dd) below, 
and the adjustment formulas are given by equations ( 51 ) . 



Sketch (dd) 
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'(<P,) = 

l x 


<Pi = 


Pboo’^o) - AC ®a’<°) “ ^ pAna’CO - ^’(0) 
<Pc 


(51) 


O&r 


Oa^. 


where 


ACn^’Co) - Cm^' (o) - Cn^^'fo) 
ACn^' (a a ) = C m c c 0 '( a a 0 ) - c ma 1 '( CI a 1 ) 


and the subscripts o and l. indicate values corresponding to the linear- 
ized theory and the adjusted theory, respectively. Notice also in 
sketch (dd) that the area corresponding to the adjusted value of 0^' 
is given by the following expression: 

= Pny^' + ^)“ a &i ACn^' (ff a ) (52) 

. The parameter Cm^j for pitching about the leading edge has been 
computed, using equations (35) to ( 52 ) for the same family of airfoils 
as that discusse'd previously and the results added to those obtained 
for Cmq’ (eq. (37))* The total damping moment Cm* + C^' is shown 
plotted as a function of Mach number in figure 7* it is evident from 
examination of figure 7 that the effects of thickness are small and are 
in the destabilizing direction with increasing thickness ratio. It 
should be noted that this result does not agree with those of Jones and 
of Wylly in references 25 and 26 which indicate that the effect of thick- 
ness is extremely stabilizing. In reference 27# however. Van Dyke takes 
issue with the results of references 25 and 26 , and offers an alternate 
solution correct to the second order in thickness. From this result, 
the following relationship for the damping-moment coefficient for pitch- 
ing about the leading edge can be deduced: 


r 


r . r, rxp[Mo 2 N- 2 ] i i 

i 

i 

i 

[ 1+ { 6 }J 


(53) 


where 


Cmq/ +Cmd, 


is the damping -moment coefficient obtained from the 


-*o 


linearized theory and 
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N = 


7+1 Mb' 

2 p 2 


X = 1 


for biconvex airfoils 


X = 3 


for double-wedge airfoils 


The trend of the damping-moment coefficient with thickness obtained from 
this result agrees closely with that obtained in the present paper. 


EXPERIMENT 
Wind Tunnel 


The experimental investigation of the damping- in-pitch character- 
istics of the various model configurations was conducted in the 
Ames 6- by 6-foot supersonic wind tunnel. In this tunnel, the Mach num- 
ber can be varied continuously within the subsonic speed limits 0.60 to 
choking, and at supersonic speeds from 1.20 to 2.00. The total pressure 
can also be changed at any time within the limits 3 pounds per square 
inch absolute to atmospheric. A complete description of the wind tunnel 
is given in reference 28 . 


Models 


The pertinent dimensions of the five wing-body configurations 
investigated in the present report are shown in figure 8. All the wings 
had sections which were symmetrical in streamwise planes and 3 percent 
thick. Wings having leading-edge sweep angles less than or equal to 45° 
(wings c, d, and e in fig. 8) had biconvex circular^arc sections in 
streamwise planes with maximum thickness at 50-percent chord. Wings 
having leading-edge sweep angles greater than 45° (wings a and b in 
fig. 8) had NACA 0003-63 sections in streamwise planes. The wings were 
identical in plan form and section shape to those investigated in the 
series of force tests reported in references 29 to 33* Strength consid- 
erations required that wings c, d, and e be built of steel. Wings a 
and b were built of aluminum. 

The models were fitted with bodies such that the distance from the 
apex of the body to the wing-body intersection was the same for each 
model. Also, the ratio of wing area to body maximum cross-sectional area 
was 17.9 for each model. Forward of the point of maximum radius, the 
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bodies were of tbin laminated wood construction and were identical in 
shape to those of references 29 to 33- The aluminum afterbodies were 
cylindrical in shape and terminated at the trailing edges of the wings 
to permit the models to be deflected with sufficient amplitude when 
mounted on the support structure. 


Model Support System 


The support system used in the present investigation was the same 
as that developed for the single-degree-of-freedom free-oscillation 
experiments of reference 10. The interested reader is referred to that 
report for a detailed description of the oscillation mechanism. For 
the present tests, two changes were made in the system which served to 
improve the accuracy of the data. The changes and the reasons for them 
are as follows: 

1. The vertical flexure pivots which provided vertical restraint 
in the tests of reference 10 were replaced by bearings after it was 
found that the heavier models of this investigation caused the pivots 
to twist laterally when the model was given its initial displacement. 

An undesirable yawing oscillation was thereby induced which destroyed 
the linearity of the system. Installation. of the bearings entirely elim- 
inated this yawing tendency. The damping due to friction in the bearings 
was somewhat larger than that of the flexure pivots, but by frequent 
checking of the wind-off damping and regular replacement of the bearings 
it was possible to maintain the tare damping well within acceptable 
limits . 


2. A stiffening strut was installed between the sting support and 
the tunnel ceiling in order to remove the possibility of coupling between 
the model oscillation and the resonant mode of the sting support. The 
influence of the strut on the aerodynamic forces at subsonic speeds was 
investigated by removing the strut and recording data for one of the 
models at a frequency sufficiently below the sting resonant frequency to 
avoid excessive sting vibrations. Results obtained with the strut 
installed and removed agreed within the experimental accuracy. 

By effectively eliminating both yawing tendencies of the model and 
vibrations of the sting support, it can be seen from conparison of the 
present results with those of reference 10 that the deviations between a 
number of observations at a given Mach number were markedly reduced. 
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Scope of Tests 


Investigation of the damping-in-pitch characteristics of the various 
model configurations was conducted over a supersonic Mach number range 
of 1.20 to 1.90 and, where possible, over a subsonic Mach number range 
of 0.60 to 0.90. The Reynolds number for the major portion of the tests 
was held constant at 1.6 x 10 s per foot. Additional data for one of the 
models was obtained at Reynolds numbers of 1.0 x 10 s per foot and 
3-2 x 10 6 per foot. 

All the tests were conducted with the models at a mean angle of 
attack of 0 °, the angle being measured from the mean line of the sting 
support to the axis of the test section. For each oscillation, an ini- 
tial displacement of ± 5 ° was imparted to the model by means of the 
pneumatically actuated pawl arrangement described in reference 10 . 

Other important variables are listed below for the five wing-body 
configurations investigated: 


Model 

Range of moment 
of inertia, 
slug - ft 2 

Axes of 
rotation, 
% M.A.C. 

Range of reduced 
frequency 

A = 2 triangular 

0 . 053 ^ - O.O 589 

35, to 

0.012 

- 0.037 

A = 3 triangular 

.0368 - .0402 

25, 35 

.018 

- .050 

A = 4 triangular 

.0491 - .0512 

25, 35, 45 

.007 

- .o4i 

A = 3 swept 

.0549 - .0576 

25, 35 

.011 

- .025 

A = 3 nnswept 

.0422 - .0488 

20 , 35 , to 

.007 

- .025 


Reduction of Data 

The technique used in this investigation for reducing an experi- 
mental oscillation-decay record to the damping-in-pitch coefficient 
Cmg + was identical to that described in reference 10. Briefly, 

this technique involves plotting the envelope of the oscillation-decay 
curve against time on semilogarithmic graph paper . If the damping of 
the system is a linear function of the angular velocity, this plot will 
be a straight line, the slope of which is proportional to the damping 
term. Calibrations for the model moment of inertia and the damping due 
to internal friction, and standard measurements for the density and 
velocity of the air stream then enable one to reduce the aerodynamic 
damping term to coefficient form. The aerodynamic restoring moment 
can be obtained from measurements of the frequency of oscillation and 
a calibration for the static spring constant of the system. Results for 
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the restoring moment about two axes of rotation a known distance apart 
then permit calculation of the lift-curve slope Cr . 

Gu 


Corrections to Data 


No wind-tunnel- wall corrections were made to the subsonic results 
for any of the measured stability derivatives. To the author's knowl- 
edge, no corrections have been developed for application to the forces 
acting on a finite-span wing oscillating in a compressible fluid. 
Tunnel-resonant-frequency effects were investigated by use of refer- 
ence and it was found that the range of frequencies used in this 
investigation was well below the lowest calculated tunnel resonant fre- 
quency. In view of the fact that the subsonic damping moments measured 
for the A = b triangular wing in the present investigation agreed 
reasonably well with the results for a similar wing obtained during an 
investigation in the Ames 12-foot wind tunnel (ref. 35) > it is believed 
that the effect of the tunnel walls on the damping in pitch was not 
significant. 

A subsonic correction could have been made to the static parameters 
and C^. However, since this correction would have been very small, 

and in view of the uncertainty involved in applying static corrections 
to the results of dynamic measurements, it was decided not to make the 
correction. 

The effects of constriction of the flow at subsonic speeds by the 
tunnel walls were taken into account by the method of reference 3 6. 

This correction amounted to, at most, a 2-percent increase in the Mach 
number and in the dynamic pressure over that determined from a calibra- 
tion of the wind tunnel without a model in place. 

For the tests at supersonic speeds, no corrections were required to 
be made on either the aerodynamic measurements or the air-stream 
properties. 

The reader is referred to reference 10 for a discussion of the 
precision of the data. Since the method of reduction of the data was 
the same as that of reference 10, the uncertainty in the recorded value 
of a given damping coefficient remained of the order of ±0.02. However, 
by virtue of the elimination of yawing tendencies in the model and sting- 
support vibrations, the standard deviation of a number of observations 
at a given Mach number was reduced from ±0.06 to ±0.03* 
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RESULTS 

Damping Coefficients 


The "basic experimental results for the damping coefficient Cniq+Cm^ 
for the five wing-body configurations investigated are shown in fig- 
ures 9 to 13« In cases where data could be obtained at subsonic as well 
as supersonic speeds, these results are shown plotted through the Mach 
number range 0.60 to 1.90 for ready comparison of the magnitudes of the 
coefficients in the two speed ranges. Also shown plotted in figures 9 
to 13 are the theoretical results for the damping coefficient Cm™ + Cm^ 
at supersonic speeds for the wings alone. In addition, theoretical 
results for the triangular wing -body combinations, considering inter- 
ference effects, are shown in figures 9 to 11. The sources of these 
theoretical results are references 11, 16 through 23, and 37- The exper- 
imental results for subsonic speeds are not compared with theory, since 
a rigorous linearized theory for triangular, swept, or unsvept wings 
oscillating in a subsonic compressible flow has not yet been developed. 

A very approximate analysis for Cm^ and Cm& for the aspect ratio h 
triangular wing was presented in reference 10 . In this development, the 
parameter was approximated using the incompressible theory of ref- 

erence 8 for the elliptic wing. However, a more rigorous analysis for 
the elliptic wing in compressible flow (ref. 11) has indicated that the 
result of reference 10 for Cm^ may be incorrect. The good agreement 
between the experimental subsonic results of reference 10 and the sub- 
sonic theory of reference 10 may be fortuitous. 


Ihe results for the three triangular wings (figs. 9 to ll) and the 
swept wing (fig. 12) show that in the supersonic speed range, in every 
case, the linearized first-order-in-frequency theory provides a reliable 
guide for obtaining both the magnitude and trend with Mach number of the 
damping coefficients. Of particular importance are the theoretical 
predictions of ranges of Mach number over which the wings could experi- 
ence negatively damped oscillations. The experimental results for 
the A = 3 triangular wing (fig. 10) are inconclusive in this regard 
since the theoretically predicted unstable range occurs In the low super- 
sonic Mach number range, 1.0 to about 1.08, where it was not possible to 
obtain data. The predictions by the theory of larger unstable ranges 
for the A = 4 triangular wing and the A = 3 swept wing were borne out 
by the experimental results, as shown in figures 11 and 12 by the exper- 
imental points plotted below the abcissas. In the absence of available 
theoretical results for unswept wings, the experimental results for the 
unswept wing are compared with theoretical damping coefficients for a 
rectangular wing of aspect ratio 3 (fig* 13 )* It is seen from examina- 
tion of figure 13 that the theoretical prediction of a rather large 
region of instability for an axis of rotation at 0.20c (fig. 13(a)) Is 
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borne out by the experimental results but that the theory fails to pre- 
dict the reversal in trend of the damping coefficients at low supersonic 
Mach numbers for axes of rotation closer to the aerodynamic center 
(figs. 13(b) and (e)). It is clear from these results that the effects 
of taper can be important and should be considered for an accurate theo- 
retical appraisal of the damping-in-pitch characteristics of unswept 
wings. 

It should be mentioned that at a Mach number of 0.90, the damping 
of the A = k triangular wing (fig. 11) and the A = 3 triangular and 
swept wings (figs. 10 and 12) were highly nonlinear with angle of attack, 
being stabilizing for angles of attack of 5° and decreasing steadily to 
zero at about 1°. The dashed lines for this Mach number shown in fig- 
ures 10, 11, and 12 are intended to indicate the range of thiB nonlinear 
variation. Further discussion of the phenomenon is withheld to a later 
section of this report. The above-mentioned difficulty was not encoun- 
tered with the A = 2 triangular wing, and the value of 
at Mq = 0*9 for this wing (fig. 9) represents the damping throughout 
the amplitude range. Subsonic results for the un swept wing could not be 
obtained due to the highly erratic behavior of the aerodynamic restoring 
moment which generally was so strongly destabilizing as to counterbalance 
the spring restoring moment and force the model against its stops. 


Reynolds Number Effects 


For the A = 2 triangular wing at supersonic speeds (fig. 9) > 
damping coefficients have been obtained at three Reynolds numbers rang- 
ing from 1.18 X 10 6 to 3*77 X 10 e , while at subsonic speeds (fig. 9) they 
have been obtained at 1.18 X 10 6 and I .89 X 10 s . In both cases, it Is 
evident that at least within the range of Reynolds numbers at which data 
were recorded there is no significant effect of scale. It is reasonable 
to assume that a similar lack of scale effects exists for the other wings 
as well. 


Aeroelaetic Effects 


In figures 14 to 23, the static parameters and Cp^ for the 

five wings are compared with force-test measurements of similar wings 
(refs. 29 to 33). Such a comparison is useful as an indication of the 
Importance of aeroelastic effects on the aerodynamic properties of the 
less rigid models of the present tests. The fact that at supersonic 
speeds there Is excellent agreement between the two experiments for the 
lift-curve Blope Cj^t whereas the present results for are 
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consistently smaller, is attributed primarily to the difference in the 
body shapes of the models of the two investigations. In the force-test 
investigation, the afterbodies were carried beyond the trailing edges of 
the wings, which would tend to shift the centers of pressure of the wing- 
body combinations farther rearward while leaving the total lift essen- 
tially unchanged. From this reasoning, it is concluded that aeroelastic 
effects on the aerodynamic coefficients of the present tests are of 
secondary importance. At subsonic speeds, there are differences in mag- 
nitude of both and Cn^. However, in view of the fact that no large 

effects of aeroelasticity were found at supersonic speeds where the 
dynamic pressure is greatest, it is believed that the differences 
in and Cn^ at subsonic speeds are not due to aeroelastic effects. 

It must therefore be stated that at the present time the differences at 
subsonic speeds are not understood. 

Transfer of Axes 


Sufficient information has been given in figures 9 to 23 to permit 
the transfer of the damping-in-pitch results at supersonic speeds to 
axes of rotation other than those used during the tests. Such a trans- 
fer, if valid, of course greatly increases the range of applicability 
of the given data. The validity of the transfer equation (given in 
appendix A) depends primarily on the invariance of the stability deriv- 
atives with changes in angle of attack and angular velocity and could be 
best checked by experimental means. Checks on the method are provided 
in figures ll(e) and 13 (c) for the A = 4 triangular wing and A = 3 
unswept wing, respectively. Shown in these figures are the experimental 
damping coefficients obtained during the investigation and coefficients 
for the same axis position computed, by means of the transfer equation, 
from the results of experiments at two other axis positions. It is 
evident that the transferred results portray the actual data with suffi- 
cient accuracy to establish the validity of the method. 


DISCUSSION 

Effect of Aspect Ratio 


It will be remembered from the discussion in the theoretical section 
of this report that the linearized theory indicates that a decrease in 
aspect ratio has a highly stabilizing effect on the damping in pitch of 
a triangular wing with subsonic leading edges. This indication is con- 
firmed in figure 24, wherein the damping in pitch is shown as a function 
of Mach number for the three triangular wings having their axes at 0.35c. 
It is noteworthy that the A = 2 wing shows no tendency toward dynamic 
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instability in the Mach number range 1.2 to '1.9* that for Mo<1.2 the 
trend of the damping coefficients of the A = 3 wing is toward insta- 
bility, and that a range of instability exists for the A = 4 wing. It 
is of further interest to note that for Mach numbers greater than 1.67, 
in which range the ..leading edges of the A = 3 and A = 4 wings are 
supersonic, the damping coefficients for these wings are essentially the 
same. This serves as a partial confirmation of the theoretical predic- 
tion of reference 19, wherein it is shown that the damping coefficients 
of triangular wings having supersonic leading edges are independent of 
aspect ratio. 

A more comprehensive way of comparing the tendencies of the wings 
toward dynamic instability would be to show the entire ranges of axis 
positions and Mach numbers over which dynamic instability is possible, 
that is, to plot the stability boundary curves for the three wings. 

In the Mach number range 1.2 to 1.9* however, only the A = 4 wing has 
a region of instability and this curve is shown, compared with the theo- 
retical curve, in figure 25. 

The preceding comparison (fig. 24) was useful primarily for verify- 
ing the theoretical prediction regarding the role played by aspect ratio 
in determining the tendency of a triangular wing toward dynamic instabil- 
ity. A comparison of the actual, magnitudes of the damping moments was 
masked, however, by the fact that the coefficients for the three wings 
were referred to their own respective characteristic lengths. In the 
following comparison, an attempt is made to overcome this difficulty by 
posing the question: Given three triangular wings of aspect ratios 2, 3 > 

and 4, of equal area, and required to have the same restoring moment 
(in ft-lb) at a Mach number of 71. 4, how do the physical magnitudes of 
the damping moments compare? For this purpose, the aspect ratio 3 wing 
is chosen as a standard, required to have a static margin of 0.05 
at M 0 = 1.4, and the damping moments of all the wings referred to the 
M.A.C. of the aspect ratio 3 wing. The comparison on this basis is 
shown In figure 26. It is evident from examination of figure 26 that on 
this basis also, the effect of decreasing the aspect ratio Is beneficial 
to the damping in pitch, the aspect ratio 2 wing having the largest 
damping moments throughout the range of test Mach numbers, and the aspect 
ratio 4 wing the smallest. 


Effect of Plan-Form Shape 


For this comparison, use Is made of the results for the aspect 
ratio 3 triangular, swept, and unswept wings. First, the prediction 
made in the section entitled "Theory" regarding the tendencies pf 
three similar wings toward dynamic instability is examined by comparing 
their stability boundaries. From the preceding section, it will be 
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recalled, however, that no region of instability exists for the A = 3 
triangular wing in the Mach number range 1.2 to 1.9* It is apparent, 
therefore, that as predicted by the theory the triangular wing is 
superior in this regard. The stability boundary curves for the swept 
and unswept wings, with axis position referred to the leading edge of 
the M.A.C., are compared in figure 27. Again, the theoretical predic- 
tion is borne out for, as seen in figure 27 , a considerably larger region 
of instability exists for the unswept wing. 

Next, the magnitudes of the damping moments are compared by choosing 
the A ** 3 triangular wing as a standard, with a static margin of 0.05 
at Mo = 1.4, and, as in the previous section, requiring that the wings 
have equal restoring moments at this Mach number. The comparison is 
shown in figure 28. On this basis, it is again apparent that the trian- 
gular wing is the superior wing, followed in order by the swept and 
unswept wings. 

Notice in figures 26 and 28 that since the damping coefficients of 
all five wings investigated are referred to the M.A.C. of the A = 3 
triangular wing and since all wings have equal restoring moments at 
Mo = lA, the results of figures 26 and 28 may be compared directly. 

Here it appears that the A <= 2 triangular wing is the most desirable 
wing from a longitudinal -dynamic-stability standpoint, and the A « 3 
unswept wing the least desirable. As mentioned previously, however, 
it should not be inferred that the result of this comparison implies 
a recommendation for the use of triangular rather than unswept wings. 
Obviously, the addition of tail surfaces could alter considerably the 
relative damping-in-pitch merits of the wings investigated. 


Effect of Thickness 


It has been shown by the results of experiments (figs. 9 to 12) 
that the supersonic linearized potential theory provides a reliable 
basis for predicting the low-frequency damping- in-pitch characteristics 
of thin wings of finite span. Considerable doubt has been shed, however, 
over the ability of the linearized theory to predict the dynamic behavior 
of wings that cannot be classified as thin. In recent years, a number 
of reports have been issued (refs. 25 , 26, and 38 ) which indicate that 
second-order thickness effects have a profoundly stabilizing effect on 
the damping in pitch of two-dimensional wings. In particular, refer- 
ence 26 indicates that by increasing the thickness of an infinite-span 
wing from zero to only 4-l/2 percent, the rather large region of insta- 
bility predicted by the linearized theory is completely el iminat ed. In 
contrast to these results, the approximate analysis of this paper 
(see Theory) indicates that thickness effects are relatively small, the 
destabilizing effect on canceled by a stabilizing 
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effect on Cm^. in view of these conflicting results, it is of some 
interest to obtain additional experimental data which would serve to 
clarify the issue. Of particular interest are the effects of thickness 
for finite-span wings. In this case, a limited comparison can be made 
between the experimental results for the 3-percent A = 4 triangular 
wing of this report and the results of reference 10 for two 6-percent 
triangular wings having the same aspect ratio. Aside from the differ- 
ences in profile and small differences in the dimensions of the bodies, 
the models and test conditions of the two experiments were almo st iden- 
tical. The comparison is shown in figure 29* It is seen that for this 
wing, the effects of thickness are small and are in the direction of 
decreasing stability with increasing thickness. Of primary importance 
is the fact that all the wings exhibit regions of Mach number over which 
dynamic instability is experienced. 

Although admittedly the experimental evidence in figure 29 is not 
conclusive, nevertheless in the light of these results it is difficult 
to accept those results which indicate that large stabilizing effects of 
thickness are present on infinite-span wings, especially In view of the 
fact that the behavior of a triangular wing with supersonic edges is in 
many ways quite similar to that of an Infinite-span wing. 


Effect of Nonllnearities 


From examination of the experimental oscillation-decay records at 
both low subsonic and all supersonic speeds. It was -evident that the 
assumption of a linear second-order system was well justified. Aside 
from a moderate amount of scatter, a decay curve could be fitted by an 
exponential curve over a wide range of amplitudes, verifying that the 
damping parameter was essentially independent of amplitude. At high sub- 
sonic speeds, however, this linear behavior was no longer true for most 
of the wings tested. As an example, figure 30 shows the oscillation- 
decay record for the A = 3 triangular wing at Mq = 0.90. It is seen 
that after an initial displacement the amplitude of oscillation quickly 
dies out to a low level, but there sustains itself indef ini tely. 

Similar phenomena have also been reported in the forced oscillation 
experiments of reference 35 for a wing-body combination having a trian- 
gular wing of aspect ratio 4 and in flight tests of various high-speed 
aircraft, particularly tailless aircraft. The deleterious effects of 
such sustained oscillations on the qualities of an aircraft as a gun 
platform or from a structural standpoint are obvious. It is therefore 
of considerable importance to gain an understanding of the phenomenon, 
with a view toward advancing means of eliminating it. 

At the present time, insufficient Information, either theoreti ca l 
or experimental, is available to enable a complete and authoritative 
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description of the physical event; a tentative hypothesis is therefore 
advanced instead, which would appear to fit the observed facts. It must 
he emphasized, however, that in the light of additional information, the 
assumptions and conclusions drawn here may require considerable revision. 

Specifically, the observed facts are these: (l) The damping 

coefficient is highly nonlinear, being negative or stabilizing for angles 
of attack greater than about 1°, and zero at 1°. The possibility that 
the nonlinearity is caused entirely by scale effects and/or wind-tunnel- 
wall interference is ruled out by the fact that similar behavior is 
observed in full-scale flight tests. (2) The frequency of oscillation 
is essentially constant throughout the amplitude range, implying that 
the restoring moment is a linear or nearly linear function of the angle 
of attack. 


From these observations, it is hypothesized that the characteristic 
differential equation governing the motion is of the form 

la + cp(a,d) a + K^a = 0 

where cp(a,a) represents the damping coefficient as a nonlinear function 
of the angle of attack and angular velocity, and I and K are constants. 
It is therefore necessary to search for a mechanism that can affect the 
damping moment to a much larger degree than it does the restoring moment. 
In this regard, it is pertinent to first review some of the character- 
istics of wings in steady transonic flow. It is well known that at high 
subsonic Mach numbers, the essential feature of the flow is the appear- 
ance of shock waves on the wing as the speed of air over the wing surface 
exceeds the local speed of sound. These shock waves cause large changes 
in the pressure and also promote flow separation near the shock wave, 
due to the presence of adverse pressure gradients. The effects of these 
disturbances on the pitching-moment characteristics of an unswept wing 
of aspect ratio 3 are shown in reference 39 and reproduced as curve (a) 
in sketch (ee) . Notice that 
the slope of the pitching 
moment is positive for smal 1 
angles of attack, indicating 
that the aircraft would be 
statically unstable in that 
range. It was also illus- 
trated in reference 39 that 
by reducing the curvature of 
the wing profile near the 
trailing edge, thereby 
reducing the tendency toward separation, the positive slope through zero 
could be eliminated. This result Is shown as curve (b) in sketch (ee) . 
By virtue of the beneficial effect of reducing separation, it is argued 
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here that in the absence of separation and strong shocks on the basic 
wing, the variation of pitching moment with angle of attack would be 
nearly linear, as shown by curve (a) in sketch (ff), and that the effects 

of these disturbances can 
be lumped into an additional 
moment variation, curve (b), 
such that the sum of (a) 
and (b) then gives the 
observed nonlinear variation. 

For triangular wings, 
the reversal in slope of the 
pitching moment at s mal l 
angles of attack does not 
appear (see refs. 29 > 30, 
and 31) • This, however, does not preclude the presence of the nonlinear 
variation due to separation. For a triangular wing with center of grav- 
ity at, say, the quarter M.A.C. point, the slope of the pitching -moment 
curve can be quite large. As seen in sketch (gg), the linear varia- 
tion (a) can then greatly outweigh the effect of the nonlinear curve (b), 
with the result that a scarcely evident nonlinearity appears in the 
combined curve. 






Sketch (gg) 


The preceding examples were taken from the 
results of static measurements, where the angle 
of attack is simply the angle between the chord 
line of the wing and the free-Btream direction. 
When the wing is oscillating, however, each 
point | measured from the axis of rotation 
experiences an additional angle of attack 
<x|/V 0 , due to the angular velocity. Then if 
the steady- state pitching -moment variation 
caused by the separation (curve (b) in 
sketch (gg)) is approximated in the range 

0<a<a max b y 


C m (a) = aa - ba s 


in the unsteady case it is. 


Cm 


+ 




( 5 *) 
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where now | is assumed to be the distance from the axis to the point 
at which the additional lift due to separation is concentrated.® Now, 
since for slow frequencies <x£/V 0 is much smaller t han a, all but 
first-order terms in a, are neglected so that equation ( 5*0 becomes 

C m = a [a-ba 2 ] + — [a-3ba 2 ] (55) 

v o 

The characteristic equation of motion then becomes 


Ici - — ^ [a-3ta 2 ] - a -^Cm^+Ca-ba ] — 0 


where [Cm^+Cm^] and Cm^ are the (constant) stability derivatives which 
would be present in the absence of separation. Thus, it appears that 
although the effect of the nonlinearity on the restoring moment can be 
insignificant, the possibility still remains that the damping moment can 
be profoundly affected in the event that [Cmq+Cm^] is sufficiently small. 
Notice that when [Cniq+Cm,^] + 2|a/c is greater than zero and if the non- 
linear term in the restoring moment can be ignored, equation ( 5 6 ) takes 
on the form 

a - 25a(l-pa 2 ) + K 2 a =0 (57) 


Equation ( 57 ) is then recognized as being the well-known Van Der Pol 
equation of nonlinear mechanics. It iB evident that for sma ll values 
of a the damping term is negative, leading to a divergent oscillation, 
whereas for larger a, the dampin g ter m is positive. A stable regime 
therefore will exist where a = y i/m- and oscillations of either large 
or small amplitude will converge to that regime. 


Thus, by the assumptions of an additional nonlinear pitching-moment 
variation caused by flow separation and a sufficiently small damping 
moment in the absence of separation, the observed phenomenon of a small - 
amplitude self-sustained oscillation can be hypothetically explained as 
being due to the destabilizing effect of the nonlinearity in the restor- 
ing moment on the damping moment. 


3 It is recognized that usually the nonlinear effects of shock-wave 
boundary-layer interaction and flow separation are markedly reduced 
during nonsteady motions. However, for the low-frequency oscillations 
of the present tests (of the order of one cycle per 100 chord lengths 
of travel) , it is believed that the effects of separation, though per- 
haps less severe than for the steady case, are nevertheless still 
present. 
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CONCLUSIONS 

The results of a theoretical and experimental investigation of the 
single-degree-of-freedom damping in pitch of a series of low-aspect-ratio 
wing-body combinations made at subsonic and supersonic speeds lead to the 
following conclusions: 

1. In the Mach number range 1.2 to 1-9* theoretical and experimen- 
tal values of the damping-in-pitch coefficient C^ + C^ were in good 
agreement for triangular wings of aspect ratios 2, 3 1 and 4, and swept 
wing of aspect ratio 3. 

2. Theoretical predictions of the existence of ranges of center- 
of -gravity positions for a range of Mach numbers greater than 1.2 over 
which dynamic instability may be expected for the aspect ratio 4 trian- 
gular wing and the aspect ratio 3 swept and unswept wings were confirmed 
by the experimental results. 

3. The prediction by the theory of the beneficial effect on the 
damping in pitch at supersonic speeds of a reduction in aspect ratio waB 
borne out by the results of experiments for the triangular wings having 
aspect ratios 2, 3> and 4. In the Mach number range 1.2 to 1-9, the 
aspect ratio 2 wing had the largest damping moments throughout the range 
of Mach numbers and the aspect ratio 4 wing, the smallest. 

4. Experimental results for the stability boundaries in the super- 
sonic speed range of three wings of aspect ratio 3 having triangular, 
swept, and unswept plan forms confirmed the theoretical prediction 
regarding the relative magnitudes of the region of Mach number and 
center-of-gravity position in which dynamic instability could be experi- 
enced. In the Mach number range 1.2 to 1.9 > no region of instability 
existed for the triangular wing. The region Of instability for the 
unswept wing was considerably larger than for the swept wing. 

5. The effects of profile thickness on the damping in pitch at 
supersonic speeds of triangular wings of aspect ratio 4 and thickness 
ratios of 3 percent and 6 percent were found to be small and In the 
direction of decreasing stability with increasing thickness ratio. 

6. The occurrence at a Mach number of 0.9 of small-amplitude self- 
sustained oscillations of the triangular wing of aspect ratio 4 and the 
triangular and swept wings of aspect ratio 3 was attributed to a desta- 
bilizing effect on the damping moment of nonlinearities in the aero- 
dynamic restoring moment. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif. 
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APPENDIX A 
TRANSFER OF AXES 


The problem to be solved may be posed as: Given the damping 

coefficient and the restoring-moment coefficient about 

two axes of rotation, a known distance apart, find the damping coefficient 
about an arbitrary axis. 


For single-degree-of-freedom rotary oscillations about axis (l) the 
damping coefficient is written as 

+ Cm^J = j^Cm^ + Cm •. + 2x 1 ^Cm a J - + c Ig.J - 2x 1 2c I tt ( A1 ) 

where the axis (o) is chosen as a reference axis, £ 0 ^ + Cm^J is known, 

and Xi is the known (nondimens ionalized) distance of axis (l) ahead of 
axis (o). For pitching about axis (2), the damping coefficient is 


where 


h + + Gn^l + 2x 2 - xJ~C L(i + C^l - 2x 2 C-^ 

J 2 L J Q L J 0 L Jo 


(A2) 


and X 2 are known. 


The relationship of the axes (l) and (2) to the reference axis (o) is 
shown in sketch (hh) . Note that if the reference axis (o) is chosen to 
be coincident with the axis passing through 
the aerodynamic center, then the quantities 


+x 


2xi 


St 


and 2x, 


Sc, 


and (A2) are zero. Since 
are known, the parameters 


in equations (Al) _ 
o 

C % 1 St 
S 




(!) (2) ( 0 ) 
Sketch (hh) 


Ct._. and may be determined from 


[' 


Sx 


Sx = 


SI " 


- l 2 


Si 


Xi - x 2 


Si 


+ X 1 C I 


-*1 


(A3) 
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With Cj^ and 
tions (Al) and 



determined, it is now possible to solve equa- 


(A2) simultaneously. 


Let 



Then 


i-b - — x 2 ] + d[xi 2 - x 2 2 ] 


Ct + Ct . 


o 


J q. JL a, 

I 

Cm q + 


x l ~ x 2 


(A4) 


a — x a [c-e] + dx x 2 = f 


The damping coefficient about an arbitrary axis, with x referred to 
axis (o) and measured as positive ahead of axis (o), is therefore, 


+ Sx 


f + xc - xe 


x 2 d 


(A5) 
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Mach number, M. 

Figure 2.- Comparison of theoretical single -degree -of -freedom 
short -period pitching stability boundaries at supersonic 
speeds for three wings of aspect ratio 3, having 
triangular, swept, and rectangular plan forms. 
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Figure J.— The effect of frequency on the theoretical single- 
degree- of - freedom short -period pitching stability 
boundaries at supersonic speeds for a triangular wing 
of aspect ratio 4 . 
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Figure 4. - The effect of frequency on the theoretical single- 
degree -of - freedom short -period pitching stability 
boundaries at supersonic speeds for a rectangular 
wing of aspect ratio 4. 
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Figure 5. - The variation with reduced frequency of the 
single - degree - of- freedom rotary damping - moment 
coefficient for the two-dimensional wing at Mach 
numbers O, 0.60, and 0.80. Axis at leading edge. 




--1 



i.u Li i.z i.o i.*t t.o i. o 

Mach number, M 0 


/./ 


Figure 6.- The effect of thickness on the damping -in -pitch coefficient 
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Figure 7. — The effect of thickness on the damping -in-pitch coefficient Cm g + C m ^. for a two- 
dimensional wing having a symmetrical biconvex parabolic -arc section. Axis at leading edge. 
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Equation of fore bodies: r * r 0 [l-(l- ffj* 

S>=2* 

/ * 50 * 

0 <x <25 * 

Equation of afterbodies : r • r„ 



(a) A* 2 triangular wing (NACA 0003-63 sections). 


Figure 8 Sketch of models investigated . 
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Mach number, M a 


(a) Axis at 0.25c. 


Figure iO.— Experimental damping -in -pitch coefficients for the wing- 
body combination having a triangular wing of aspect ratio 3. 
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(b) Axis at 0.35c. 


Figure tO.— Concluded. 
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Mach number , M e 
(a) Axis at 0.25c. 

Figure II.— Experimental damping- in- pitch coefficients for the wing - 
body combination having a triangular wing of aspect ratio 4 . 
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— Continued . 
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Mach number, M 0 
(a) Axis at 0.25c. 


Figure !2.- Experimental damping-in- pitch coefficients for the wing- 
body combination having a swept wing of aspect ratio 3 . 
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Mach number, Af 0 
(a) Axis of 0.205. 


Figure 13.— Experimental damping -in- pitch coefficients for the wing - 
body combination having an unswept wing of aspect ratio 3. 
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Figure !3. — Continued . 
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Figure / 4.— Experimental pitching -moment -curve slopes for the wing- 
body combination having a triangular wing of aspect ratio 2 . 
Axes of rotation at 0.35c and 0.45 c . 
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Figure 16.— Experimental pitching -moment -curve slopes for the w/ng- 
body combination having a triangular wing of aspect ratio 4. 
Axes of rotation at 0.25c, 0.35c, and 0.45c . 
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Figure IQ.— Experimental pitching -moment -curve slopes far the wing- 
body combination having an unswept wing of aspect ratio 3. 
Axes of rotation at 0.20c, 0.35c , and 0.40c, 
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Mach number , M„ 


Figure 19.— Experimental lift -curve slopes for the wing -body 
combination having a triangular wing of aspect ratio 2. 
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Figure 22.— Experimental lift -curve slopes for the wing- body 
combination having a swept wing of aspect ratio 3. 
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Figure 23.- Experimental lift- curve slopes for the wing -body 
combination having an unswept Wing of aspect ratio 3. 
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Figure 25.— Comparison of theoretical and experimental single- 
degree- of -freedom short -period pitching stability 
boundaries for a wing- body combination having a 
triangular wing of aspect ratio 4 . 
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Figure 26. — Comparison of magnitudes of experimental rotary damping- moments for three 

wing - body combinations having triangular wings of aspect ratios 2, 3, and 4. 
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Figure 27. — Comparison of experimental single - degree - of - 
freedom short - period pitching stability boundaries for two 
wing — body combinations having swept and unswept wings 
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Figure 28. — Comparison of magnitudes of experimental rotary damping moments for three 
wing -body combinations having triangular, swept, and unswept wings of aspect ratio 3. 
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Figure 29. - Comparison of experimental damping -in -pitch coefficients for three wing - body 
combinations haring triangular wings of aspect ratio 4 and thickness ratios of 3 and 6 percent. 
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Figure .30. - Oscfflograph record at Mach number 0.90 for the wing -body combination 
having a triangular wing of aspect ratio 3. Axis at 0.256. 






